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1. INTRODUCTION 



Supersymmetric gauge theories are interesting theoretical laboratories. They are rich 
enough to exhibit most of the quantum field theory phenomena, yet they are rigid enough 
to contain a lot of exactly calculable information [1] [2] [3] [4] . They embed easily into string 
theory, and provide an exciting arena for the search of string/gauge dualities [5] [6] [7] [8] . 

In the past year a lot of progress has been made in understanding some of this rich 
structure using direct field theoretic techniques ([9] [10] and refs therein), in the case of the 
theories with extended supersymmetry, and in [8] in the case of M = 2 susy broken down 
to M = 1 by the superpotential. 

In particular, a connection between four dimensional gauge theories and two dimen- 
sional conformal field theories seemed to appear. Some earlier indications for such a 
connection were observed in the study of M = 4 super- Yang-Mills theory [11] [12]. 

In this paper we shall make this connection more transparent in the case of M = 2 
theories. 

We shall also derive, by purely field theoretic means, via direct instanton calculus, the 
solution for the low-energy effective theories, proposed by Seiberg and Witten in 1994 [13] 
and further generalized in [4] [14] . 

1.1. Notations 

Throughout the paper we denote the colour indices by the lower case Latin letters I, m, n = 
1, . . . , AT, the vevs of the Higgs field by 



will occur often. 

The gauge group G = U(N) : the group of gauge transformations on R 4 which extend 
smoothly to S 4 will be denoted by Q. Its normal subgroup Q^, which consists of the gauge 
transformations, trivial at oo, will play a special role. 




a = — 




the dual vevs by 



£ = diag(fi,...,fjv). 



The vector 
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1.2. Organization of the paper 

The section 2 is addressed at physicists, who want quickly jump onto the subject. 
It contains previously unpublished details on the noncommutative regularizations of the 
theories we consider, as well as systematic introduction into fi-backgrounds. In this paper 
we shall not touch upon the recently revived [15] C-backgrounds of [16] [17], postponing 
the explanations of the relations between O and C to some future work. 

Section 3 starts setting up the stage for the mathematical problem and the reviews 
the ingredients needed for its solution. We present the formula for the partition function 
of the pure Af = 2 gauge theory in the O-background, as a sum over random partitions. 
We also begin to formulate the same sum in terms of random paths, which arise as the 
boundaries of the Young diagrams, representing the random partitions. Mathematically 
oriented reader can skip the section 2 and proceed directly to 3 (we also recommend [18] [19] 
for orientation). 

Section 4 attacks the problem of the calculation of the prepotential. Physically it has 
to do with the limit, where the O-background approaches the flat space. In the random 
path representation, this limit is the quasiclassical limit, which can be evaluated using the 
saddle point method. This is exactly the idea of our derivation. We extensively discuss the 
equations on the minimizing path, and their solution. We find that the solution is most 
simply described in terms of the Seiberg-Witten curve. 

Section 5 explains the fermionic representation of the partition function in the special 
O-background, preserving twice as much supersymmetries compared to the generic O- 
background. Even though the formalism of free fermions is well-known to mathematicians, 
under the name of the infinite wedge representations of gl(oo) algebra, we supply the 
necessary details. We find that a certain transform of the partition function can be written 
as a matrix element (current conformal block) of the exponentials in the U(N) currents 
on the sphere S 2 . 

Section 6 begins our quest for generalizations. We discuss the softly broken M = 4 
theory, i.e. M = 2 gauge theory with the matter hypermultiplet in the adjoint represen- 
tation. As most of the steps are similar to the pure gauge theory case we move faster, 
and write down the expression for the partition function in terms of partitions, paths, and 
chiral fermions, which this time live on elliptic curve, determined by the microscopic gauge 
coupling. Again we perform the saddle point evaluation, and find that the prepotential is 
encoded, as conjectured by Donagi and Witten [20], in the spectral curves of the elliptic 
Calogero-Moser integrable system. 
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Section 7 considers gauge theories with matter in the fundamental representation, and 
the gauge theories with the tower of Kaluza-Klein states, coming from the compactification 
of the five dimensional theory on a circle. 

Section 8 presents our conclusions and the discussion of the unsolved problems. 

1.3. Acknowledgments 

We thank A. Braverman, R. Kenyon, A. Losev, G. Moore, A. Morozov and S. Smirnov 
for useful discussions. Research was partly supported by NSF under the grant DMS- 
0100593 (AO), by the Packard Foundation (AO), by P$$H under the grant 01-01-00549 
(NN), by the grant 00-15-96557 (NN) for the scientific schools, and by Clay Mathematical 
Institute (NN). 

AO thanks Institut Henri Poincare for hospitality. NN thanks Institute for Advanced 
Study, New High Energy Theory Center at Rutgers University, Kavli Institute for Theoret- 
ical Physics at University of California Santa Barbara, Forschungsinstitut fur Mathematik, 
ETH Zurich, University of Sao Paulo, and ICTP, Trieste, for hospitality during preparation 
of the manuscript. 

Note added: H. Nakajima informed us that he and K. Yoshioka found another proof [21] 
of our main theorem of the section 4, which relates the partition function, prepotential 
and the Seiberg-Witten curves. They are using the blowup techniques. 
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2. Af = 2 GAUGE THEORY, DEFORMATIONS AND BACKGROUNDS 

In this section we remind the construction of [9] [10] and also provide quite a few new 
details. We consider Af = 2 supersymmetric gauge theory with the gauge group U(N). 
We study Euclidean path integral with the fixed vev of the adjoint Higgs field, on R 4 
with certain nonminimal couplings. Most of the interesting dynamics in such theories is 
associated with instanton [22] effects [23]. 

2.1. Lagrangian, fields, couplings 

The simplest way to write down the action of the four dimensional super- Yang-Mills 
theory with extended supersymmetry is to use dimensional reduction of the higher dimen- 
sional minimal supersymmetric theory. 

In particular, starting with the six dimensional Af = 1 SYM we arrive at: 




Here, Ar,0, <fi are the components of the six dimensional gauge field, which decompose 
as the four dimensional gauge field and an adjoint complex Higgs field (or two real Higgs 
fields); A a i are Af = 2 gluions - the pair of the four dimensional Weyl spinors, transforming 
in the adjoint representation of the gauge group. 

The indices et,f3 = 1,2 correspond to the doublet of the SU(2) L , a, (3 = 1,2 is that 
for SU(2)r, while i, j = 1, 2 are the internal indices, which reflect the SU{2)j R-symmetry 
of the theory. They are raised and lowered using the SU(2) invariant tensor £12 = —£21 = 
Space-time Lorentz indices will be denoted throughout the paper by 
the upper case Latin letters /, J, .. = 1, 2, 3, 4. The Pauli tensors, relating the spinor and 
vector indices, are d; a( i. 

In (2.1) we have put the bare coupling constant go and the bare theta angle -$o- The 
bare coupling corresponds to some high energy cut-off scale fx. 

The action (2.1) is the limit of the six dimensional action of the theory put on a six 
manifold T 2 x R 4 with the standard flat product metric. 
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2.2. fi-background 

However, in going from six to four dimensions one may have started with a nontrivial 
six dimensional metric. In particular, by reducing on the two torus one may have considered 
R 4 bundles with nontrivial flat SO(4) connections, such as the space Mq with the metric: 

ds 2 = Adzdz + giJ (dx 1 + V 1 dz + V* d£} (dx J + V J dz + V J dz^ (2.2) 

with V 1 = QjX J , V 1 = OjX J , and the area A of the torus to be sent to zero. For [O, O] = 
the metric (2.2) is flat. 

However, for O ^ the background (2.2) will break all supersymmetries. Indeed, 
the spinors e^,eai generating supersymmetries of the (2.1) are the components of the six 
dimensional Weyl spinor. In order to generate the symmetry of the theory on a curved 
background the spinor must be covariantly constant. In our case this means that the 
spinor, viewed as four dimensional Weyl spinor, should be invariant under the holonomies 
around the two cycles of the two-torus, which is possible only for discrete choices of O, O. 

2.2.1. fi-background in the physical formalism 

Fortunately there exists a continuous deformation of the M = 2 theory preserving 
some fermionic symmetry (which also deforms along the way). The trick is to use the 
R-symmetry, which is manifest in the four dimensional theory — the SU(2) group, which 
acts on the internal index i = 1, 2 of the gluinos A. 

Namely, in addition to the nontrivial metric (2.2) we turn on a Wilson loop in the 
R-symmetry group, which should compensate some part of the metric induced ho lo no my 
on the spinors. 

As a result, in the limit A — > with O, O fixed the action (2.1) gets extra terms: 

AL = tfjL { p J + UjL ( P J + OSof Lf K JL (2.3) 

where 

4 1} J = Jd 4 xVG (x J {TvFtkD^ + ^TtX^DjX^ + G JK a^Tr[A di , A?]) (2.4) 

where the tensor (7j K = — = o 3 ^ = \sikjl^ JL is the 't Hooft projector: so(4) = 
su L (2) © su R (2) - su R (2) = ((I) ® {\)) sym . 
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and, for future use: 

13 1 / pa f3a\ 

Finally, 

Lf K 3L = Jd 4 xVG x J x L G MN TrF IM F KN (2.5) 



2.2.2. O-background in the twisted formalism 

Supersymmetric gauge theories with extended supersymmetry can be formulated in 
a way which guarantees the existence of a nilpotent symmetry in any curved background. 
This formulation, sometimes called twisted, or topological, or cohomological, makes use of 
the R-symmetry of the theory. The coupling to the curved metric is accompanied by the 
coupling to the R-symmetry gauge field, which is taken to be equal to the corresponding 
projection of the spin connection [24] [25]. In this way, the fermions of the pure gauge theory 
become a one-form ip, a self-dual two-form x, and a scalar rj. The bosons Aj and are 
not sensitive to the twisting. 

The advantage of such a formulation of the theory is the clear geometric meaning 
of all the terms in the action. It is well-known, that M = 2 super- Yang-Mills in the 
twisted formulation provides an integral representation for the Donaldson invariants of 
four-manifolds (for G = SO (3)). In our study the gauge theory lives on R 4 which is 
boring topologically. However, we study something different from the Donaldson theory, 
due to Q-background. It corresponds to the K-equivariant version of Donaldson invariants, 
K = Spin{&) being the group of rotations. 

In general, if the metric G = Gijdx I dx J of the Euclidean space-time manifold has 
some isometries, one can deform the standard Donaldson- Witten [24] [26] action, by cou- 
pling it to the isometry vector fields V and V, which should commute [V, V] = 0. Explicitly, 
the action of the theory is given by: 

T ■-— (-|TrF * F + Tr {D A <f> - t v F) * (D A - tyF) + |Tr[0, 0] 2 vol 5 ) 



2^ 

+Tr (x {D A ^) + + VD A ^ + x * L V X + *L V V + V> * If!*) 

+Tr (x * [(f), x] + V * [</>, v]+ip* [0, 

+— TrF A F 
2n 



(2.6) 



On R 4 , we take, as above, 



V 1 = tfjx- 7 , 



v 1 = nW 



(2.7) 



with 





f ° 


ei 
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n IJ = 


-ei 
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0/ 



(2.8) 



2.3. On supersymmetry 

For completeness, we list here the formulae for the super symmetries of the actions 
(2.1)(2.3)(2.6). 

2.3.1. Supersymmetry of the physical theory 

The supercharges transform in the representation (2, 1, 2) © (1, 2, 2) of the Lorentz x 
R-symmetry group SU(2)lxSU(2)rxSU(2)i. Introducing the corresponding infinitesimal 
parameters: £a and £ d the supersymmetry transformations can be written 



SX ai = ai^Fjj + iC a D + iV^a^DjBe^d 



'a/3 
U 



(2.9) 



SB = V2( ai \ ai 
SB = V2(t^ 

with the auxiliary field D = [B, B]. 

2.3.2. Twisted superalgebra 

The twisted formulation of the theory is achieved by replacing the Lorentz group 
K = SU(2) L x SU(2) R G K x 5?7(2)j by another subgroup of K x SU(2)j, namely 
SU{2) L x SU(2) d , with SU(2) d being diagonally embedded into SU(2) R x SU{2)l In 
other words, the internal index i is identified with another SU(2)r index, a. The fields 
are redefined, according to 



, x a/3 -T-a/3 T af3 



(2.10) 
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Similarly, the supersymmetry parameters C a become C 7 , Cdi become ( IJ and (. Of course, 
( IJ is self-dual. 

The supersymmetry algebra becomes: 



SAi = - KiJip 3 - i(iV 

Hi = +C J Ffj + KiD + i(D I( j) 
Sxu = (F+ - KijD + i (C [7J Dj] 
S v = ( IJ F 

8<j> = CV/ 



(2.11) 

5 V = ( IJ F IJ -i(D + ( I D I <j ) 



5<P = (V + C U XiJ ■ 

The geometric meaning of the transformations 5 is the following. The space of the 
fields Ai, 4>, ipi, represents the ^-equi variant de Rham complex of the space A of gauge fields 
on R 4 , together with the ingredients needed to construct a Mathai-Quillen representative 
of the Euler class of a certain infinite dimensional bundle over A, and the projection 
form[27] associated with the projection A — > A/Q. 

The space R 4 is hyperkahler, and possesses an action of the group R 4 of translations. 

The transformation generated by ( is the ^-equivariant de Rham differential. The 
transformations, generated by ( IJ correspond to the ^-equivariant dx,j,K differentials, 
corresponding to the three complex structures on A induced from the complex structures 
on R 4 . Finally, Q 1 correspond to the operators lq; of contraction with the vector fields on 
A, induced by the vector fields 

d 



dx 1 

generating translations of R 4 . 

2.3.3. Supersymmetry of the fi-background 

The transformations, which generate the symmetry of the O-background utilize the 
rotational symmetries of R 4 . To the vector field 

v - ^ 

there is an associated vector field on A and the associated operation of contraction on the 
£-equi variant de Rham complex. In terms of the transformations S of twisted theory these 
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are simply the transformation (2.11) with the space-dependent transformation parameter 
The theory (2.3) in invariant under the transformation, generated by: 

(CX I X IJ ) = (CV I (x)C0) 

The supercharge, generating this transformation, will be denoted as Q. 
2.4. Noncommutative deformation 

The theory (2.1) allows yet another deformation, which we shall implicitly use to simplify 
our calculations. Let Q IJ = —<3) JI be a constant Poisson tensor on R 4 . The theory on R 4 
can be deformed to that on the noncommutative space, R@. A naive way to define this 
deformation is to replace all the products of functions (or components of various tensor 
fields) in (2.1) by the so-called Moyal product: 



f*g(x) =exp-0 iJ 



f(x + Og(x + v) (2-12) 

77=C=0 



2 dt 1 dr] J 

More conceptual definition goes as follows (cf. [28] [29]). Consider the theory (2.1) di- 
mensionally reduced to zero space-time dimensions. We get some sort of supersymmetric 
matrix model. Replace the matrices by the operators in the Hilbert space 7i. Explicitly, 
we get the theory of six bosonic operators 

X 1 , $,$; 7=1,..., 4 (2.13) 
and eight fermionic ones (we use twisted formulation): 

* J , ??, x n T = \s ijkl xkl\ 7,J=1,...,4 (2.14) 

where we lower indices using some Euclidean metric gij: 

Xkl = 9ki9ljX IJ 7 

The action reads: 



L = 



4 g 



^Tr H [g IK g JL ([X 1 , X J ] [X K , X L ] + X IJ [X K , V L ] + X /J [$, X KL \) 
+ 2g u (ryLY 7 , * j ] + [¥, X 1 ] [$, X J ] + V 1 [¥, * J ]) 
+ ([$,$] 2 + r ? [$,r ? ])j 

+ Tt H ( iBjj [X I ,X J ] + ^- e IJKL [X 1 , X J ][X K , X L ] + L 
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(2.15) 



The action (2.15) leads to the equations of motion which do not depend on the parameters 
entering the last line in (2.15), i.e. Bu 7 i} , L . 



A special class of extrema of the action (2.15) with all the fermionic fields set to zero is 
achieved on the operators X 1 , $, $ obeying: 

[X 1 , X J ] = iQ IJ I, Q IJ = 2g 2 g IK g JL B KL 

(2.16) 

[X 1 ,®] =0= = [$,$] 

Let us fix a standard set of operators x J in H obeying 

[x 7 ,x J ] =i@ IJ 1 (2.17) 
The algebra (2.17) may be represented reducibly in H. In general, for nondegenerate /J , 

H = H ®W (2.18) 

where Ho is irreducible representation of (2.17), and W is a multiplicity space, which we 
shall assume to be a finite dimensional Hermitian vector space, of complex dimension N. 
From now on, x J will denote the operators on Ho- The corresponding operators on H will 
be denoted as x J <g> 1 N . Expanding: 

X 1 = x / <g> 1 N + z6 /J Aj(x), $ = 0(x) = 1 ® a + 0oo(x) , * J = /J ^j(x) (2.19) 

where 4>oo{x) — > 0, x — > oo, we arrive at the naive formulation, with the metric on R 4 , 
given by: 

G /J = 9 klO IK Q JL 

$o=$ Pf (6) (2 .20) 



^ = ^pf(e)v^ 



9o g 2 

In (2.19), both Aj = Aj\ ln <g> E ln and <j> are valued in the N x N matrices, anti-Hermitian 
operators in W. To arrive at the ordinary gauge theory, we should take the limit 0^0, 
while keeping Gu, g%, $q finite. The curvature F of the gauge field A is given by: 

[x 1 , x J ] = i<?> IJ + e IK e JL F KL (2.2i) 
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2.4.1. Combining and 

The universal gauge theory (2.15) can also be subject to the nontrivial O-background: 

^ TTn(^9iK9jL^ [X 7 ,X J ][X K ,X L ]+ 

+ x IJ [x K , * L ] + X IJ ([*, x KL ] + ^lx ML ) 

+ 2g u X 1 ] + U[x L ^j ([$, X J ] + n J K X K ^j + 

+ rilX 1 , V J ] + # J * J ] + fi^ L ) 



+Tr n [ iB IJ [X I ,X J ] + ^ e IJKL [X I ,X J ][X K ,X L ]+L 



(■ 

(2.22) 

The vacua of the theory (2.22) are given by the operators which solve slightly different 
equations then (2.16). Namely, the condition on $ is now: 

[X 1 , $] = QjX J , [X 1 ,®] = q!jX j (2.23) 

which is consistent with [X 1 , X J ] = z6 7 J only under certain conditions on 6, O. Namely, 
assuming nondegeneracy of 0, with u> = _1 , the matrices 

£ u = Q*fu>Ki, Si j = UjUki, (2.24) 

must be symmetric. Then the solution to (2.16) is given by: 

X 1 = x 7 , $ = ±£,jx 7 x J , $ = ±£ 7J x 7 x J (2.25) 

Now, the gauge field Aj and the Higgs field </> are introduced via: 

X 1 = x 7 <g> 1 N + z0 /J Aj(x), $ = ±£ 7J X 7 X J <g> ljy + <£(x) (2.26) 
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2.4.2. Supersymmetry of the noncommutative theory in background 

We shall only write down the supercharge of our interest: 
QX 1 = # 7 , Q^ 1 = X 1 ] + n T jX J 
Qd> = , Q X IJ = H IJ 
QH" = [*, X IJ ] ~ {SlW* - ^W K ) + 
Q§ = rj, Qrj = [<£>,<&] 

Here H IJ is an auxiliary field, which is equal to [X 1 \X J ] + on-shell. 

2.4.3. Observables in the supersymmetric gauge theory: = 

In the ordinary Af = 2 supersymmetric gauge theory a special class of observables 
play an important role. They are distinguished by the property, that a certain super- 
charge annihilates them. Of particular interest for application to Donaldson theory are 
the observables which are constructed out of the invariant polynomials in the Higgs field <p 
by means of the descend procedure. In terms of the twisted fields, these observables look 
as follows. Let P{4>) be any G-invariant polynomial on the Lie algebra g = Lie(G) of the 
gauge group G. Then: 

of(x) = P($(x)) 

f BP eft P 



Of(X) = / 



d 2 P 1 <9 4 P 

i . F a AF b + _ + n . f A^Af Af 



l x 2 d(j) a d(t) h " ' 24 d(t> a d(f) b d(t) c d4> d 

Here, x, C, E, . . . , X represent a 0, 1, 2, ..,4-cycles in the space-time manifold, respectively. 
The main idea behind (2.28) is to use the fact that the supercharge Q acts on P($(x) + . . .) 
as the de Rham differential. To compare the observables of the ordinary gauge theory to 
those of the noncommutative gauge theory we shall utilize the generating function(form): 

P = P ((f>(x) + ij I {x)dx I + \F IJ (x)dx I A dx J ) G fi* (spacetime) (2.29) 

Its main property is: 

d x (O P (x)...) = 
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The noncommutative gauge symmetry does not allow, naively, for the local observables 
like in the first line of (2.28). However, it has as many gauge invariant observables as does 
the ordinary theory. For the gauge group U (N) the invariant polynomials can be expressed 
as polynomials in the single trace operators 

P n (0) = Tr0 n (2.30) 

It is convenient to introduce the character: 

Pp(4>) = Tre** = V { -^fp n (2.31) 

The analogue of (2.28) (2.29) is the following closed form on R 4 : 

P(3 (x, dx) = pe-P" J dHd^K e-Mw'+Oidx 1 ) ^jm*,*) ( 2 .32) 

where 

= $ + K/X 7 + tf 7 * 7 + ±tf / tf J [X 7 ,X J ], u = \uijdx I Adx J (2.33) 
The closedness of Op^ is proved with the help of the symmetry acting on k: 

5$ = k, 5k = (2.34) 



2.4.4. Observables: O^O 

The chiral observables get deformed when O ^ 0. First of all, in the ordinary gauge 
theory, the generating form Op becomes equivariantly closed under correlator: 

(d + i V )0 P = (2.35) 

In the noncommutative gauge theory we get the same statement, but now we need to 
modify the symmetry 5 to its equivariant analogue: 

5$! = kj, Sk t = -ttjilj (2.36) 
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and change (2.32) to: 

a=l,2 Cq ^ 

(2.37) 

= \ (uJudx 1 A — £[jx I x J ) 

In what follows we consider either O-observable, or integrated 4-observable. In the first case, 
we set x = (to be at the fixed point), and get the integral over the k of Tr- W e l/3 *^' K ' ) . 
In the second case, we get the integral over the k of Tr- H e lf3 ®(' & ' K ' with an extra 5- 
invariant Gaussian factor. Actually, (2.36) implies that the integral over k is localized 
at k = # = 0, so that both 0- and 4-observables are equal, up to a e a -dependent factor: 

° { pU = n ^ - 4 ^* 

a=l,2 

R4 a=l,2 Ca 

The observables (2.38) are most natural in the five dimensional gauge theory compactified 
on the circle of the circumference f3. The periodicity of Op p in $ has there a simple origin 
- large gauge transformations [14] . We shall return to this theory in the section 7. 



L 
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3. GAUGE THEORY PARTITION FUNCTION 



In the ordinary gauge theory the natural object of study is the Euclidean path integral 
over the field configurations, such that 

<j)(x) — > a, x — > oo (3-1) 

We are aiming to calculate the path integral in the fi-background, with the boundary 
conditions (3.1) (and, in fact, using noncommutative deformation): 

Z(a,ei,e 2 ;A) = J DADtJjDi]DxD$D$ e~ -fc 4 L v ^ d * a! (3.2) 

In the noncommutative gauge theory the condition (3.1) is phrased differently. Think of 
$ as of the element of the Lie algebra of the group Go ~ U(H) of gauge transformations 
(note that normally one takes Qq = PU(H)). The identification H = C N ® Ho allows 
to view <E> as an N x N matrix of operators in Ho- The condition (3.1) is now stated as 
follows: 

$ = a ® l Wo + ± 1 N ® ^ / jx / x J + p(x) (3.3) 
where <f(x) — > as x — > oo, faster then any power of x. 



The integral (3.2) needs an UV cutoff \i. The bare couplings go and $o are renormalized, 
thus generating an effective scale A, 

OM OAT -^+27rii? 

A 2N _ jj 2N ( , IT ( 3 4 ^ 



Remark: On ?7(1) factor. The gauge theory with the gauge group U(N) has an interesting 
pattern of the renormalization group Row. The ordinary gauge theory has a decoupled U(l) factor, 
which does not exhibit any renormalization of its coupling, and interacting SU(N) part, with the 
famous phenomenon of asymptotic freedom, reflected in (3.4). When the theory is deformed hy 
Q, the perturhative loop calculations are altered, in a Q-dependent way, as it introduces, among 
other things, an energy scale. In particular, the U(l) factor coupling constant starts running, in 
the energy range 

^ >>A>> ^ 
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where the formula (3.4) with N = 1 holds (this result is a simple generalization of [30], which 
can also be justified using [1][31]). The [/(l)^ -1 C SU(N) gauge couplings experience different 
renormalization group Bow, as they are affected by the loops of the charged W-bosons. fn principle, 
we should introduce two distinct low-energy scales, one for SU(N) part of the gauge group, another 
for U (1): A, Ao- As we shall eventually remove noncommutativity it makes sense to keep Ao ~ fj,. 
We shall encounter the ambiguity related to the £7(1) factor in the analysis of the Seiberg-Witten 
curves, fn the paper we concentrate on the SU(N) dynamics, and mostly set a = 0. 



Remark: The theory (2.1) has an anomalous U(l) g h symmetry (gh for ghost). Under this 
symmetry, the adjoint Higgs held has charge +2, <& has charge —2, the fermions ip have charge 
+ 1, XiV h& ve charge —1. In the background of the instanton of charge k the path integral measure 
transforms under the ghost U(l) with the charge —AkN. The O deformation breaks U(l) 9 h- ft 
can be restored by assigning Q, O the charges +2, —2, respectively. 



3.1. Partition function as a sum over partitions 

The partition function Z(a, ei, e 2 ; A) can be explicitly evaluated as a sum over instan- 
tons. Moreover, the fi-background lifts instanton moduli, leaving only a finite number 
of isolated points on the appropriately compactified instanton moduli space as a full set 
of supersymmetric minima of the action. The evaluation of (3.2) is then reduced to the 
calculation of the ratios of the bosonic and fermionic determinants near each critical point. 
These points are labeled by the colored partitions. 

Consequently, (3.2) is given by the sum over A^-tuples of partitions [9] [10] (see appendix 
B for the notations and definitions related to partitions). Explicitly: 

Z(a; ei, e 2 , A) = Z^a; A, Cl , e 2 ) £ A 2Ar l% c (a; e u e 2 ) (3.5) 

k 
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where: 

7 , ■, TT a l ~ a n + ei(i - 1) + 
Zg(a;ei,e 2 ) = I| ~ 

l,n;i,j a l~ a n + Z\ (i ~ k n j - 1) + t 2 {k U - j) 

tt ai - a n + ei(-z) + e 2 (j - 1) 

l,n;i,j + tiihj - i) + e 2 {j - k ni ~ 1) 

= 1 tt r {hi - knj + v{j + + bi n ) T (u(J - i) + b ln ) 

e 2Ar|£| n -A? ^ r " fc »J + ^ " + 6 '») r ^ " Z + !) + 6 J») ' 

(3.6) 

&«»=^^, 1/=-^ (3.7) 

and 

Z^a; Cl , e 2 , A) = exp I 5> £l ,e 2 (a* - a n ; A) J (3.8) 

\l,n / 

where the function 7 £l , e2 is defined in the appendix A 

Remark: The product over i,j in (3.6) in is infinite and needs precise definition. Here it 
is: Fix a pair (l,n). Consider all the factors in corresponding to l,n. Split the set of indices 
into four groups: 

z+ = s++ u s + - u S- + u S— 

S++ = l<i< £(ki), l<j< k nl } 

S+- = l<i< £(ki), k nl < j} (3.9) 

s-+ = {(i,j)\e(ki)<i, i<j<k nl } 
s— = {(i,j)\ e(ki)<i, k nl <j} 

The set S++ is finite, the set S-- contributes 1. The set S-\ contributes: 



oo .. . ... 

y in TT TT ai n + - I) + e 2 {-j) _ , 

Z + k ll 11 a ,„+ e .i(i-l)+ e .o(k,<-i)- [ ' 



ain + ei (i - 1) + £2 {ku - j) 
=i j=fe„i+i 



nn 



ain + ei (i - 1) + e 2 (&h - k n i - j) 



and similarly for S |_. See also [32]. 
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The special case 62 = —e\ = ft, i.e. v = 1 deserves a special attention. In this case, the 
expression for the partition function simplifies to: 

Z(a;ft,A) = ^A 2 ^l^(a;ft) 
k 

Z c (a;ft) = ^(a;ft)^(a,ft) 

,,2/„ *\ _ TT I Cll — Cln + h(kl,i — k n j + j — i) \ 

^ ' " < JL, V — ^CFO J ( ' ' 

i^ ert (a; ft) = exp j ^ 7R (a, - a n ; A) 

\ l,n 

where the function 7^ (x; A) is defined in the appendix A. 

3.2. Plancherel measure on partitions 

In the case when N = 1 and v = 1, the sum in (3.5) is over a single partition k and 
the weight /U^(a, ft) reduces to 

^>=n(*^^)=n^. (3-) 

j<i v 7 nek v 7 

where second the product is over all squares □ in the diagram of the partition k and h(n) 
denotes the corresponding hook-length. 

The weight /u(k) 2 is known as the Plancherel measure on partitions because of the 
relation 

M(k) = ^, (3.13) 

where is the irreducible representation of the symmetric group corresponding to the 
partition k. Just like in the classical Plancherel theorem, the Fourier transform on the 
symmetric group an isometry of the L 2 -spaces with respect to the Haar (i.e. counting) 
measure on the group and Plancherel measure on the set of its irreducible representations. 
Observe that from (3.13) it follows that 



E ^( k ) = ^ ( 3 - 14 ) 



n 

k \=n 
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The Plancherel measure is the most fundamental and natural measure on the set of 
partitions. In many aspects, the set of partitions equipped with the Plancherel measure 
is the proper discretization of the Gaussian Unitary Ensemble (GUE) of random matrices 
[33]. In particular, the integrable structures of random matrices are preserved and, in fact, 
become more natural and transparent for the Plancherel measure, see [34] [35] . Similarly, 
Plancherel random partitions play a central role in the Gromov-Witten theory of target 
curves, extending the role played by random matrices in the case of the point target, see 
[36]. 

A pedestrian explanation of the relation between Plancherel measure and GUE can be 
obtained by rewriting the weight (3.12) as the product of the Vandermonde determinant in 
the variables ki — i and a multinomial coefficients, which, of course, is the discrete analog 
of the Gaussian weight. 

As we will see below in Section 5, for N > 1 and v = 1 the partition function is again 
related to the Plancherel measure, but now periodically weighted with period N. Finally, 
the case v ^ 1 leads to the Jack polynomial analog of the Plancherel measure, see [37] . 
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4. PREPOTENTIAL 

In this section we study the limit ei,€2 — > of the partition function (3.2). In this 
limit, according to the field theory arguments [9] the partition function must behave as: 

Z(a;ei,e 2 ,A) = exp ( — JF(a; e u e 2 , A) J (4.1) 

V e l e 2 / 

where T is analytic in ei,e2 for — > 0. We prove the conjecture of [9], which identifies 
jF (a, A) = jF(a; 0, 0, A) with the Seiberg-Witten prepotential of the low-energy effective 
theory. 

4.1. Quasiclassical limit and SW curve 

By setting v = 1 we deduce from the field theory prediction (4.1) that for % — > the 
partition function (3.11) has an asymptotic expansion: 

Z(a; ft, A) = exp ^ ft 23 " 2 JF^a, A) j (4.2) 

where 

^o(a,A) = -I^(a z -a n ) 2 (log(^ T ^)-^+ f;A 2 ^ /fc (a) (4.3) 

l,n ^ ^ ' ' k=\ 

is the so-called prepotential of the low-energy effective theory. The coefficients fk for 
k = 1,2 were computed directly from instanton calculus in [38] [39], for < 5 in [9]. 

The prepotential was identified in [3] [4] [40] with the prepotential of the periodic Toda 
lattice. The arguments for this identification were indirect. In [41] the coefficients fk for 
k < 5 were computed for the Toda prepotential, in agreement with the later results of [9]. 

The purpose of this paper is to give a proof of this relation to all orders in the instanton 
expansion. 

More precisely, we shall show that: 



The set ^ai, . . . , a^; ^f^, • • • > f^") c °i nc ides with the set of periods of the differen- 
tial 



dS = z — (4.4) 



on the curve C u defined by 



A N ( w + I ) = p N ( z ), p N ( z ) = Z N + u 1 z N ~ 1 + u 2 z N ~ 2 + ... + u N (4.5) 
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Remark: Normally the factor 2ni is included in the definition of the prepotential [18]. We chose 
not to do this to get real expressions for JF for real a; 's and A. 

For generic u = (ui, U2, ■ ■ ■ , un) the curve (4.5) is a smooth hyperelliptic curve of 
genus N — 1. The differential dS has poles at two points oo± over x = oo. The homology 
group 

Hi(C u - {oo + ,oo_};Z) 

is isomorphic to Z © Z 2 ^ -1 ), the first summand being generated by the small loop 
around oo + (the loop around oo_ is in the same homology class). The intersection pairing 
identifies the dual space with 

H X (C W {oo + , oo_}; Z) , 
the relative cycle, dual to being the path connecting oo_, oo + . 
• 

Note that in the iV-parametric family of curves (4.5) there is a set of singular curves. 
The identification between the periods and (a, i ^ a ) is up to the action of the monodromy 
group. We need to fix this identification in some open domain of the parameter space. 



Consider the region Uoo of the parameter space, where u s 3> A s . In this domain the 
curve (4.5) can be approximately described as follows. Let cki, . . . , be the zeroes of the 
polynomial Pn(x): 

N 

P N (x) = ]l(x-ai) (4.6) 
l=i 

For small A for each / we can unambiguously find af 1 ~ ai, such that 

P N (af) = ±2A N (4.7) 
These are the branch points of the two-fold covering p : C u — > P 1 , P 1 3 x. 
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Let be the 1-cycles on C u which are the lifts of the cycles on the x-plane, which surround 
the cuts going from a~[ to af . Let b/ be the relative 1-cycle, represented by the path on 
C u which starts at oo_, goes to a+, and then goes on the second sheet to oo + . The cycles 
a^, and b; have canonical intersection pairing: 

a ; n b m = 5 lm . (4.8) 

The position of the a-cycles and b-cycles on the curve C u is illustrated on the Fig.l 



Fig.l The curve C u and the cycles on it, the closed ones are aj's, the noncompact ones are b;'s 

The cycles cr^ and l°° are related to a; , b m via: 

t7 0O = ^a I , r° = ^b, (4.9) 

l i 

The periods ai and ^ L are defined via: 



o, = d> dS, ^ = 2ttz / dS (4.10) 



where b/(//) is the regularized contour of integration, which connects w\ z — ll and w~ x \ 
instead of oo_|_ and oo_. 



The divergent (with fj) part of the periods is easy to calculate: 

1 dJ- 

— = 2Nn — uilog (fx) + finite part (4-H) 

2ni oai 

At the same time: 

dS = a = ^ai = -u x (4.12) 

l 

As we said above, in the sequel we set a = = u\. Also, to avoid worrying about the 
cut-off we work with the absolute cycles bz — bz+i. 
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4.2. Partitions and their profiles 



The standard geometric object associated to a partition is its diagram. In this pa- 
per, we will draw partition diagrams in what is sometimes referred to as the Russian 
form (as opposed to the traditionally competing French and English traditions of drawing 
partitions). For example, the following figure 




Fig. 2 HnarpaMMa K)Hra 



shows the diagram of the partition (8, 6, 5, 3, 2, 2, 1, 1). The upper boundary of the diagram 
of k is a graph of a piecewise-linear function fk(x) which we will call the profile of the 
partition k. Explicitly, 



/ k (x) = \x\ + 



i=i 



\x — ki + i — 11 — \x — ki + i\ + \x + i\ — \x + i — 11 



(4.13) 



The profile is plotted in bold in Fig. 2 . 

For general ei > > ei, it is convenient to extend the definition of /k(x) by scaling 
the two axes by — e\ and 62, respectively. For example, for (€1,62) = (—1, |), the scaled 
diagram of the same partition (8, 6, 5, 3, 2, 2, 1, 1) and the corresponding profile /k(^|ei, £2) 
will look as follows 
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Fig. 3 Squeezed Young diagram 



We have 



h(x\e 1 ,62) = \x\ + 



i=i 



|x + ei — e 2 ki — e\i\ — \x — e 2 &i — e\i\ 
— \x + ei — eii| + |x — eii| 
|x + e 2 - eikj - e 2 j\ — \x - eikj - e 2 j\ 

-\x + e 2 - e 2 j\ + \x- e 2 j\ 

By construction, the profile of a partition satisfies 

fL(x\e u e2) = ±1 , 
/k(a?|ei,e 2 ) > |x| , 
/k(^ki, e 2 ) = |^|, for \x\ ^> . 
We also define the profile of a charged partition 

/ a; k(^|ei, e 2 ) = /k(z - a\ei, e 2 ) 
The charge and the size are easily recovered from f a ;k(x): 

a=W dx x f". M {x\e u e 2 ) = - \l dx f' a . M (x\e u e 2 ), 
Jr Jr 



k i _ ' f j™ _2 f// , _ 1 / .,.2 
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Here and in what follows, we denote by 

4 g(x) dx = lim / g(x) dx 

JD L ^°°^°JDn[-L,L]\ S ing s (g) 

the principal value integral over a domain D C R, where sing^g) denotes the 5- 
neighborhood of the singularities of g(x). 



For a colored partition k and a vector a we define: 



N 



/ a; k(z|ei,e 2 ) = fa|ei, £2) 



1=1 

For example, for e 2 = — e\ = h, a\ = — a 2 = llh and the partition 

{(7,4, 3, 3, 2,1), (8,7,4, 4, 3,1)} 
the corresponding profile looks as follows: 




(4.17) 



-20 -10 10 20 

Fig.4 The profile of the colored partition {(7, 4, 3, 3, 2, 1), (8, 7, 4, 4, 3, 1)} 

4.3. Vacuum expectation values and resolvents 

The profile function is natural from the gauge theory point of view: the vacuum 
expectation values of the single trace operators Tr(p n have a simple expression in terms of 
the profiles: 

(Tr^)a = -=r—± tt £ A 2JV|C '^( a ; ^ A ) n[k] (4.18) 

Z(a;ei,e 2 , A) ^ 



27 



where 



On[k] = | I dxx n f'^(x\e u e 2 ) (4.19) 



The second derivative /" of a partition profile is a compactly supported distribution on R 
which for random partitions plays to the role similar the role of the spectral measure for 
random matrices. In particular, it is convenient to introduce the resolvent R(z\ei, €2) of a 

— * 

colored partition k by 

f f"Ax\e 1 ,€2) 

R{z\e u e 2 ) = \ dx — . (4.20) 



R 



X 



Note that: 

R{z\e 1 ,e 2 ) = — + 4 + 0(-^), z — ► 00 (4.21) 

In our context, the limit shapes / of partitions will be convex piece- wise analytic functions 
for which the second derivative f"(x) will be positive and compactly supported. 

4.4. Real vs. complex 

In the gauge theory the vacuum expectation values ai of the Higgs field are, in general, 
complex. In this case the profile function / a £ does not make sense. However, the resolvent 
(4.20) of a colored partition is well-defined and can be effectively used in the analysis below. 
The formula (4.19) is replaced by: 

O n [k) = l-j z "R( z \e u e 2 )dz (4.22) 

where the contour goes around z = 00. 

Complex values of ai can be reached from the real values by analytic continuation. 
By the same token, it is enough to analyze the problem for values of ai in any open set 
of R^. In particular, it is enough to analyze the problem in the asymptotic domain Woo, 
which is what we will do next. 

4.5. The thermodynamic limit 

Our strategy in extracting the ei, €2 ^0 limit of the sum (3.5) is the following. The 
typical size of the partition k contributing to the sum is of the order |k| ~ and is 
so large that the sum over the partitions can be approximated by an integral over the 
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space of continuous Young diagrams f(x). A continuous Young diagram, by definition, is 
a function satisfying the following weakening of the condition (4.15) 



f(x) = \x\, \x\ > , 
\f(x)-f(x)\<\x-y\, 

-fdxf'(x) = 0, ( 423 ) 

/ dx (f(x) — \x\) < oo 
Jn. 

Note that the condition f'(x) = ±1 is replaced by the weaker Lipschitz condition. 

We will show that this path integral is dominated by a unique saddle point (in fact, 
a strict maximum). In a more mathematical language, this means the following. By 
associating to every partition its profile, we get from (3.5) a measure on the space of 
Lipschitz functions. As ei, e 2 — > 0, these measures concentrate around a single point, that 
is, converge to the delta measure at a single function. This function is the limit shape of 
our random partition. We will construct this limit shape explicitly and show that it has a 
very simple and direct relation to the Seiberg-Witten geometry. 



We begin by observing that the Plancherel measure //(k) can be written in terms of 
profile fk(x\h) as follows 

Mk)=exp {~j tt(x\h)f£(y\h) lh (x-y;A) dxdy^. (4.24) 

More generally, we have 

Z £ (a;ei,e 2 ,A) = exp J dxdy f'^(x\e 1 ,e 2 ) ^^^,62) j Sl ,e 2 (x - y,A)^ , (4.25) 

as can be easily checked with the help of the main difference equation. Denoting the 
right-hand side of (4.25) by Zf(ei, £2, A), we have 

Z(a;ei,e 2 ,A)= £ Z f (e u e 2 ,A), (4.26) 

£ £ P 'di s cr et e 

where the summation is over the set Y^ screte of paths of the form f = f a £- 
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When ei,e 2 — > with v fixed, the size of a typical partition k; in (4.26) grows like 
|k;| ~ j^-. In this limit, the sum (4.26) looks like an integral, over the space T a of paths 
/ of the form 

JV 

f(x) = J2fi(*-ai) (4-27) 
i=i 

with each // satisfying (4.23). For this integral, eie 2 plays the role of the Planck constant. 
Our strategy, therefore, is to find a saddle point (in fact, the minimum) of the action 

*a(/) = \J dxdy f"(x)f"(y) (x - yf (log - |) (4.28) 

on the space T a of paths (4.27). The action (4.28) is the leading term as ei, e 2 — > of the 
action in (4.25): 

Z f ( a; ei, e 2 , A) ~ exp (—E A (f)) (4.29) 
Integrating by parts we rewrite (4.28) as 

W) = ~\J < (N + f\x))(N- f\y)) log (^j^jdxdy (4.30) 

which, for N = 1, reproduces the result of Logan-Schepp-Kerov-Vershik[42] [43] [44] [45]. 
Thus, we have: 

^ (a,A) = -Crit /6 r a £ A (f) (4-31) 

4.5.1. Profiles vs. eigenvalue densities 

In [9] an expression for Zfc(a, ei, e 2 ) = |k|=fc was gi ven i n terms of some contour 
integral over k eigenvalues 0j, I = 1, . . . , k. This expression follows straightforwardly from 
the ADHM construction of the moduli space of instantons, and, therefore, easily generalizes 
to the case of SO, Sp gauge groups. It is, therefore, quite remarkable that the formula 
(4.29) can be obtained directly from the contour integral expression, avoiding the actual 
evaluation of the integral (which is, of course, needed to get correclty the J^'s with g > 0). 



i=i 



Cl + ^2 



.27reie 2 P{(f>i)P{(f>i + e 1 + e 2 ) 



n 

i±j 



VIJWIJ + ei + £2j 
(<PlJ + Cl)(0JJ + £2) 



(4.32) 
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where 



P(x) = Y[(x 



a>i, 



(4.33) 

4>u = 4>i -4>j 

We now multiply Zf~ by A 2kN and sum over k = 0, 1, . . ., to get Z inst . Take the limit 
ei,e 2 — > 0. The typical k which will contribute most to the sum will be of the order 



k ~ So we introduce the density of eigenvalues: 



p{x) = eie 2 ^S(x - 



(4.34) 



i=i 



which is normalized in a /c-independent way, which nevertheless guarantees its finiteness 
in the limit we are taking. Now, the difference with the ordinary 't Hooft-like limit of 
the ordinary matrix integrals is the presence of the equal number of the 4>u terms in the 
numerator and the denominator of the measure (4.32). This will change qualitatively the 
density dependence of the effective potential on the eigenvalues, and the resulting equilib- 
rium distribution of the eigenvalues. In particular in our limit the former superymmetric 
matrix integral (4.32) scales as exp/cF as opposed to the 't Hooft's exp/c 2 F. 

Nevertheless we have a sharp peak in the measure, which justifies the application of 
the saddle point method. Indeed, by expanding in ei,e2 we map the measure in (4.32) 
onto: 



where 



A 2fcA %(a, ei ,e 2 ) 

p(x)p(y) 



exp 



ei£2 



-EaH 



dxdy 

'x^y (x-yV 
Up to the perturbative piece 



2 / dx p(a;)log 



P(x) 



(4.35) 



(4.36) 



Ln 



(ai - a n ) 2 log 



a i 



a r 



A 



the energy (4.36) coincides with the action (4.28) if we identify: 




(4.37) 

The same method applies to other theories considered in [9] [46] . 

Thus we have learned yet another interpretation of the limiting profile of the colored 
partition. In the rest of the paper we shall work with limit profiles, as the equations on 
the extremum are identical to those following from (4.36). 
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4.5.2. Surface tension 

Given a function f(x) as in (4.27), how can we extract az's? This is easy to do in the 
region where 

ai < ai +1 , 

in which case the supports of the functions f at ,u.i do not intersect. Introduce parameters 
£1, . . . , £jv which will play the role of dual variables to the charges a;. From (4.16) we have 
the following 

Proposition 1 For m <C a t+1 we have 

y]6oi = -\j cr(f'(x))dx, 
i J * 

where cr(y) is a concave, piecewise-linear function on [-N, N] such that 

o-'(y) = Ci , ye [-iV + 2(Z-i),-iV + 2Z] , 

and 

cr(-iV) = -a(iV) = -J]6- 
An example of the graph of a is plotted on the Fig. 5 




Fig. 5 Surface tension for £ = (4, 1, —2, —3) 

In Section 5, we define a dual partition function Z D A). We will see that the 
dual partition function Z D can be interpreted in terms of a periodically weighted (with 
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period N) Plancherel measure on partitions. The periodic weights in this formalism are 
precisely and the function a becomes the corresponding surface tension function. 

This story has many parallels and direct connections to the theory of periodically 
weighted planar dimers developed in [47]. 

It is a general principle that singularities of surface tension (usually referred to as 
"cusps" ) correspond to flat regions ( "crystal facets" ) of the corresponding action minimiz- 
ing shapes. In our case, the singularities of the function a correspond to gaps between the 
supports of the functions — a/) for the minimizer f*(x). 

4.6. The equations for the limiting shape 

In order to minimize £(/) with fixed a/'s, we introduce Lagrange multipliers £/, which 
we order as 

£1 > • • • > Cat , (4.38) 
and look for the maximizer of the total action 

S A (f) = SkU) + \ J *(/'(*)) dx (4.39) 

for fixed values of the £j's. To get back the minimum of £(/) with fixed ai we shall later 
perform the Legendre transform with respect to £/'s. Note that the surface tension term 
\ I a (f) i n (4.39) under the conditions a = and (4.38) can be made arbitrarily large by 
making the individual profiles fi sufficiently separated from each other. The price one pays 
for this is the decrease in the "energy" term — £a(/)> thus there is a competition between 
the two terms in (4.39) which leads to the global maximum. 

The action (4.39) is a concave functional. In fact, the first term in (4.39) is strictly 
concave (with proper boundary conditions) , which can be seen by rewriting it as a certain 
Sobolev norm, see [44]. Therefore, any critical point of the action S(f) is automatically a 
global minimizer. It is, therefore, enough to look at the first variation of (4.39). Because 
of the singularities of a, this first variation will involve the one-sided derivatives. 



Taking the first variation and integrating once by parts, we find the following equation: 

y-x 



4 dy (y-x) | log 

J y^x 



A 
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1 f"{y) = o'{f{x)) . (4.40) 



This equation should be satisfied for any point x, for which f'(x) is a point of continuity 
of a'. When 

f'{x) e {-N + 21 | I = 1,...,N- 1} 
then considering the left and right derivatives separately we obtain the inequalities 



Xf(x)e(*'(f'(x)-0),*'(f'(x) + 0)) , 



(4.41) 



where, by definition, 



[X/] (x) = 4 dy (y - x) I log 

Jy^x \ 



y-x 



A 



- 1 f"(v) 



The transform X/ is closely related to the standard Hilbert transform 



[H<7] (x) = --f 



dy 



g{y) 



Indeed, 



[X/]" = 7TH(/") 



(4.42) 



(4.43) 



(4.44) 



With = +oo, £n+i = —oo, the conditions (4.40) and (4.41) can be recast in the following 
form: 



Proposition 2 A function /*(x) is a critical point of <S(/) iff 

X/*(x)=£i whenever - N + 21 - 2 < fj{x) < -N + 21 
& > X/*(o;) > 6+i whenever //(x) = -N + 21, l = 0,...,N 



(4.45) 



In other words, the function 



<p{x) = f ir '{x) + ±[Xf 1r ]'{x) 



(4.46) 



defines a map from R to the boundary of the domain A depicted on Fig. 6. We now 
describe the solution to (4.45) in great detail. 

34 



4.7. Construction of the maximizer 

Let <&(z) be the conformal map from the upper half-plane to the domain A which is 
the half-strip 

A = {w | \%t{m)\ < N,5s(w) > 0} (4.47) 

with vertical slits along 

= -N + 2l,S{w) G [0, 77/]}, Z = l,...,iV-l. (4.48) 




Fig. 6 Conformal map for N = 3 



The positive reals rji, I = 1 . . . N — 1, are the parameters of function $. 
We normalize $ by the condition that it maps infinity to infinity and 



$( z )=iV+— log- +0(-) 

TTl Z Z 



OO . 



(4.49) 



This fixes $ up to an overall shift x ^ x + const. This ambiguity is related to the overall 
shift of the limit shape (and hence, to the overall charge of our colored partition) and is 
immaterial. 



4.7.1. Construction of the conformal map 

The map $ can be, of course, found using the Schwarz-Christoffel formula, but it is 
easier to construct is as the following sequence of elementary conformal maps. 
We take Pn{z) to be a monic real polynomial 

P N (z) = z N + ... 



35 



such that all roots of of the equation 

P N (z) 2 - 4A 2N = Y[(z - a+)(z - af) (4.50) 



N 



1=1 



are real. Let w be the smaller root of the equation 



A N [w + — ) =P N (z). (4.51) 



The function 



w 



w h-> A N ( w + — 

w 



is known as the Zhukowski function and it maps the open disk \w\ < 1 to the exterior 
of the segment [-2A iV , 2A Ar ]. It also maps reals to reals, and therefore the smaller root 
of (4.51) maps the upper half-plane to the disk \w\ < 1 with slits along the real axis. It 
remains to take logarithm to obtain the map $. Concretely, 

<&(z) = —log (w) + N , (4.52) 

where w is the smaller root of the equation (4.51) and we take the branch of the logarithm 
satisfying 

31og (w) — > , z — > +oo . 



Since 



A N 



we get the normalization condition (4.49). 

4.7.2. Example: N = 2 

As an example, consider the case N = 2 
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-2A N 2A N 




Fig. 7 A sequence of maps for N = 2 
A real quadratic polynomial 

P 2 (z) =z 2 + ... 

maps the upper half-plane to the entire complex plane with a cut along the ray 



minlVaO, +oo 

xER J 

In particular, if the min^i^f^) < —2 A 2 , that is, if both roots of 

P 2 (x) = -2 A 2 

are real then the segment [— 2A 2 , 2A 2 ] is contained in this cut. 

4.7.3. Gaps and bands 

For general N, the only difference with the N = 2 case is that the map Pn winds the 
upper half-plane N/2 times around the complex plane. Correspondingly, the root of (4.51) 
lives on a certain cover of the disk \w\ < 1, which then gets unfolded by the logarithm. 

Observe that the conformal map <&(z) has a well-defined extension to the boundary 
R of the upper half-plane 

ip(x) = $(x + zO) , 
37 



where the notation is consistent with (4.46). 

We label the roots of (4.50) so that the union of the intervals [c^ - , af] is the preimage 
of the base of the half-strip A under the map <p, that is 

N 

( p- 1 ([-N,N]) = \J[aT,a+}. (4.53) 
l=l 

These intervals are plotted in bold in the above figures. We will call these intervals the 
bands of $ and the complementary intervals — the gaps of $. 

Proposition 3 We claim that for a choice of £i's in (4.38) we can find the corresponding 
values of 77; , I = 1, . . . , N, such that the function /* satisfying 



(4.54) 

is the maximizer of the action S(f). The bands and gaps of $ will correspond to the curves and 
flat parts of the limit shape /*, respectively. They will also correspond to the bands and gaps in 
the spectrum of the Lax operator for the periodic Toda chain, as will be explained in the next 
section. 



Let us apply the Schwarz reflection principle to any vertical part of the boundary of A and 
the corresponding piece of the real axis in the domain of $. Taking the resulting function 
modulo 2, we obtain an AT-fold covering map 

N 

$ mod 2 : C \ (J[ a t , Oii~] — > {Qw > 0}/ mod 2 (4.55) 
1=1 

from the complex plane minus the bands to the half-infinite cylinder, shown schematically 
on Fig. 8. 

Omod 2 



Fig. 8 Covering half-cylinder 
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Note that this map has square-root branching precisely over the points irji, I = 1, . . . , N— 1. 
It follows that &(z) extends to an analytic function in the complement of the bands. By 
the reflection principle, the value of &'(z) on the other side of the cut is precisely —$>'(z). 
Also note that &(z) is real (resp. purely imaginary) on bands (resp. gaps). We conclude 
that 

= w W{x)_ dx = _^ R {z) (456) 

iri J R x — z iii 
where Rf+{z) is (ei, e 2 — > limit of) the resolvent of the limit profile defined by (4.54). 
We can, therefore, identify: 

Z'th band = [af, ctf] = supp/ ; (x — a/), I = 1, . . . , N (4-57) 



4.7.4. Periods 

It follows from (4.56) that 



1 



ai = - — : <p zRf i< (z)dz = <p dS 



(4.58) 



where we used (4.52), and, in agreement with (4.4), we have defined Seiberg-Witten dif- 
ferential: 



1 dw 
db = z — 

Z7TZ W 



(4.59) 



Taking the average of &(z) on two sides of the cut in (4.56), we get 



H%'(x) = -$f(p'(x) . 



(4.60) 



Using (4.44), we conclude from (4.60) that 



[X/*]'(x) = -n^ip(x) + const . 



We claim that this constant is, in fact, zero. Indeed, we have 

PW(s) = - f log f"(y) dy = 2iVlog + O Q } , x^oo (4.61) 
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Comparing this to (4.49), we conclude that 



[X/*]'(a;) = -icQKp(x) . (4.62) 

Since, clearly, ^sip(x) > and ^sip(x) vanishes on the bands, it follows from (4.62) that 
X/* is monotone decreasing and constant on the bands. It is also clear that the function 
(4.54) is monotone increasing and constant on the gaps, where it takes the values 

-N + 21, l = 0,...,N. 

This verifies the conditions (4.45) with £/'s being the values of X/* on the bands. 

To recover the £/'s we just integrate [X/*]' along the gap, where 

i$$(p(x) = <fi(x) + const . 

We find 



6+i - & = - 




integrating by parts and using the vanishing of ^s<p(x) at the endpoints of a gap. In terms 
of the Seiberg-Witten differential and the cycles b; the last relation reads as: 

£/ - 6+i = 2m I dS (4.63) 

Jbi— b; + i 

This together with the overall constraint ^2£i = fixes the values of the &'s. 

4.7.5. Completeness 

We now show that by a suitable choice of the slit-lengths r\i we can achieve any value 
of the parameters (4.38). The period map 

R^o" 1 3 (77i, . . . , 77iv-i) » (£i > • • • > £iv) , (4.64) 

is an continuous (in fact, analytic) map of open sets of R^ -1 . Because of the uniqueness 
of the maximizer this map is one-to-one. Hence, if we can additionally show that it 
maps boundary to boundary, it will follow that this map is onto. 
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It is clear that if 77^ — > for some / then £1 — — > 0. Suppose that for some 
n=l,...,N, 

i] n = msixtji — > +00 

It is clear from the electrostatic interpretation of the conformal map $ that the tip of the 
n-th slit will not be screened by other slits, that is, for some constants 81, 82 > we have 



/ 



dx > 82 

, ,a~],Q<S>(x)>(i ln -S 1 ) 



whence 

fj - = / dx > 8 2 (r] n - 81) -> 00. 



4.7.6. Periods and the prepotential 

We now perform a check. We consider the partials 

gg(A) 

and relate them to the A-periods a; (4.58) of the differential dS. 
We have 



(4.65) 



(A 



(4.66) 



because any infinitesimal change in /* can only decrease the value of <S(/*) which forces 
the variation of 5 due to the change in /* to vanish (this is essentially a usual argument 
about the variation of the critical value of the action with respect to the parameters of 
action). The rest is trivial: 

SS(f+) =\ J [8a](W) dx=-\J [8a]'{fJ) xff(x)dx = 

x f*"( x ) dx = (4.67) 

l Ja i 

which is exactly what we wanted, given (4.52) (4.59). 
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Thus, 



^o(a, A) = 


-&(/*) 


7~~> / \ 7 

RfAzjdz = 


71 / \ 

oflog (w) 


ai = <f> dS 




J SLl 




& = 2wi f 


dS 



(4.68) 



The integral in the last line is to be understood with the cut-off. Again, in the SU(N) 
theory, where only the differences ai — am, and — £ m make sense, the cut-off never show 
up. It is the Cheshir cat smile of the noncommutative regularization. 

4.8. Lax operator 

The Seiberg-Witten curves arising from conformal maps $ can be parameterized as 
the spectral curves in the periodic Toda chain corresponding to the real initial conditions. 

Consider the infinite Toda chain with particle coordinates qi and momenta pi. Make 
it periodic by imposing the constraints 



q i+N = qi - N log A 



(4.69) 



for all i. The Lax operator of this periodic Toda chain is a discrete Schrodinger operator 
of the form 



/ Pi 



,91-92 



L(w) 



wA N e qN ~ qi \ 



,91-92 



P2 e 



,92-93 



92-93 



,9-iV-l — 9JV 



\ w - 1 A N e qN - qi 



,9iv-i — 9jv 



(4.70) 



Pn / 



where w is Bloch-Floquet multiplier. 

The integrals of motions are summarized by the spectral curve, which is the curve 
defined by the characteristic polynomial 



det(z - L(w)) = P(z) — A N I w + 



w 
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Here P(z) is a monic polynomial of degree N. Observe that all roots of the polynomials 

P(z) ±2A N = det(z-Lfrl)) (4.71) 

are real because L(±l) is a real symmetric matrix. This fact plays an important role in 
the dynamics of the periodic Toda lattice. 

The band and gaps of the map $ are precisely the bands and gaps in the spectrum of 
the associated periodic discrete Schrodinger operator L on Z 

[Lf] (<) = e «-i-*/(i - 1) + Pi f(i) + e*-*+i/(i + 1) , i e Z . 

Indeed, by (4.69) L commutes with translation operator T 

[Tf] (i) = f(i + N) 

and (4.70) is the restriction of L onto the iV-dimensional w-eigenspace of T. In a band, 
the Bloch-Floquet multiplier w is complex number of absolute value 1 and hence L has a 
bounded eigenfunction. 

It can be shown that all curves of the form (4.51) with iV real ovals arise in this way. 
This is similar to the result of [47] that the spectral curves of periodically weighted planar 
dimers parameterize Harnack (also known as maximal) plane curves. Incidentally, they are 
also M-curves, not just because they are used in the M-theory construction of the gauge 
theory [7] but also, and mostly 1 because they have exactly real N ovals. 

For example, taking the lattice at rest leads to Chebyshev polynomials. In this case 
all gaps shrink to points. 

4.9. An SU(3) example 

Here is an example of a limit shape /*. Take 

P(z) = z 3 -4z. 

To visualize the curve w + — = P(z), let's look at the plots of 3?(it?) and 5s (w) for z G [—3, 3] 
plotted in the Fig. 8 in bold and normal, respectively 

1 We should apologize to M-theorists for the fact, that M-manifolds were introduced long time 
before the lid sugra was invented 
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-3 
Fig.9 




4.10. Higher Casimirs 

Just as in [10] we could deform the theory by adding arbitrary higher Casimirs to the 
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microscopic prepotential: 



^ M 2 + ^r,n(jTr$ 



n J=l 



(4.72) 



The deformations by the single trace operators are especially simple, as they would lead 
to the modification of the action S(f) by the purely surface term: 

S(f;r n ) = -£ A (f) + - J dx a(f') + - J dx /"(*) E rfc fq^ ( 4 - 73 ) 



Note that T\ shifts log (A). Presumably the critical point of (4.73) would be a solution of 
some generalization of Whitham equations [48] [40] [18] [49] [50] [51]. 
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5. DUAL PARTITION FUNCTION AND CHIRAL FERMIONS 



In this section we show that a certain transform of the partition function (3.2) has 
a natural fermionic representation. The physical origin of these fermions is still not com- 
pletely clear. One way to understand them is to invoke the Chern-Simons/closed string 
duality of Gopakumar-Vafa [52], and then the fermionic representation of the large N 
topological gauge theories [53] . Another possible origin is through the M-theory fivebrane 
realization of the gauge theory. Mathematically this story is very much related to the 
Heisenberg algebra representation in the cohomology of the moduli space of the torsion 
free sheaves on CP 2 , studied in [54]. 

5.1. Dual partition function 

Let £i, . . . , £jv be the complex parameters, 

l 

Consider 

z D {i- P - h,A) = 

Z{h{ Pl + Pl );h,A) expf ^5>'6) ^ 

Clearly, 

Z D (Z;p + N;h,A) = Z D (Z;p;h,A) (5.2) 

(shift all p/'s by 1). Thus, there are essentially N partition functions one could consider. 
They are labeled by the level 1 integrable highest weights of sl/v- Moreover: 

z D (Z;p + i;h,A) = z D (t; P ;h,A) (5.3) 

where 

£ + = (6,6,---,6v,£i) (5-4) 

To extract from Z D the partition function of interest we perform a contour integral. In the 
search for prepotential we are actually interested in the extremely high frequency Fourier 
modes of Z D , as we want JF (a; A) as a function of finite ai = h(pi +pi), with % — > 0. This 
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means that the inverse Fourier transform can be evaluated using the saddle point, which 
we already analyzed. 
Clearly, 

oo 

Z D (^p;h,A) =exp^^- 2 ^;p,A) (5.5) 

3=0 

where T§ is in fact p-independent, and is given by the Legendre transform of JF : 

^ D (^;A) = z^6az+^o(a;A), 6 = ^ > ( 5 - 6 ) 
i.e. must be given by the b; periods of the differential eiS 1 . 

Remark: Note that in this section ditfer by the factor i from the £j 's of tiie previous section. 



5.2. Free fermions 

Introduce A" free chiral fermions ip^: 

l 



r6Z+2 



1 



r6Z+2 



which can also be packed into a single chiral fermion $ 

{* r ,$ s } =5 r+s 

•w = v (f ) * , » w = E (t) * ^ 

in the standard fashion [55] : 

y N (r +Pl )=i>r\ ^N{r- Pl )=¥r ] (5-9) 
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The operators \I/ r , \& s act in the standard fermionic Fock space H (sometimes called an 
infinite wedge representation). It splits as a sum of Fock subspaces with fixed U(l) charge, 
defined below: 

H = ® pe zH[p] 

Introduce affine U(N) 1 currents, which act within H\p] for any p: 

J'»(z) =: V^ (n) :=E^T E (5-10) 

Here we normal order with respect to the vacuum |0), which is annihilated by 

tf r |0) = 0, r>0 
* s |0) = 0, s>0 

which is equivalent to 



V? } |0> = 0, r>0, 

$ z) |0>=0, s>0 
The normal ordered product is simply: 



# r * s , s > 
- # s # r , r>0 



WrWs > - ^ m V?\ ^>o 



(5.11) 



(5.12) 



(5.13) 



One can also introduce vacua with different overall £7(1) charges: 

* r |p) = 0, r>p, 

$ a |p) = 0, s > -p (5.14) 

b) e W[p] 

It is also useful to work with \I> and the corresponding C/(l) 1 currents: 

J(z) :(z) = -^J2 J lm (z~ N )z l ~ m (5-15) 

l,m 

and Virasoro generators: 

L = r : ^r^-r ■ (5.16) 

Note: 

^ob) = ^-b> (5-17) 
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5.2.1. Bosonization 

It is sometimes convenient to work with the chiral boson: 

<f>(z) =q-ij log (z)+i V -J n z~ n 

n 

dcf) =zd z (j) = —ij 

[q, M = i 

xfj( z ) = : e **(*) : , q( z ) = : e -**(*) : 
We will also use a truncated boson, with the zero mode q removed: 



1 



(z) = -ij log (z) + i Y] ~Jn z 



(5.18) 



(5.19) 



The space H\p] is actually an irreducible representation of the Heisenberg algebra generated 
by J: 



H\p] = Span ni nk>0 J- ni ...J- nk \p), Jo 



= P 



H[p] 



We also recall: 



* : n\p] -> n\p + 1], * : w[p] -> n\p - 1] 



(5.20) 



(5.21) 



5.3. Dual partition function as a current correlator 

We claim: 



Z D (Z;p;h,A) = (p| e */' ftB +W-'W e /* ,&ff W-'WA 2I '«» e -*/ ,ftB -W J W|p> (5.22) 



where the matrices if are given by: 



N 



E+(z) = zE N ' 1 + J2 El ~ 1,l i 



1=2 



H(z)=J2*iiE 



i,i 



N 



E_(z)=z- 1 E 1 ' N + J2E 1 ' 1 ' 1 



1=2 
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(5.23) 



5.4. Affine algebras and arbitrary gauge groups 

Note that (5.22) can be written using the Chevalley generators e^, fa, hi of the afhne 
Lie algebra sljy: 

Z D ^;p;?i,A) = {u n ,A 2ho e^u h ) v , 



uh = exp ^ w ( 5 - 24 ) 

AT-l 

;=i 

where V UJp is the integrable highest weight module with the highest weight vector vq (an- 
nihilated by the simple roots e{) and the highest weight u p , p = 1, . . . , N. 
The formula (5.24) has an obvious generalization to any simple Lie algebra g. Conjec- 
turally, it gives the partition function in the O-background in the M = 2 gauge theory 
with the gauge group, S'-dual to G (Langlands dual). 

5.5. Dual partition function and gl(oo) 

We now return to the Ajy-i case. 
In the language of the single fermion \1/ the formula (5.22) reads as: 

Z D (^p;h,A) = (p\e^e^A 2L °e^\p) (5.25) 

where is a diagonal matrix: 



H £ = -Vf i : # r #_ r : 

^ % Z^ s (r+5)modjV r r 

r 

Clearly, [H^,L ] = 0. The formula (5.25) expresses the dual partition function as an 
average with the Plancherel measure of the iV-periodic weight e H « . 

Remark: On Toda equation 

The function Z D obeys Toda equation (cf. [56]): 

< g(A) io g (z>» = ^fr+gr 1 ' (5.26) 

In fact, the formula (5.25) identifies Z D as the tau-function of the Toda lattice hierarchy, 
thanks to the results of [57] [36] [58], with the specific parameterization of the times. We 
hope to return to this property of the partition function in a future publication. 
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5.6. The U(l) case 

In order to understand (5.25) (5.22) we consider first the case N = 1. As the Plancherel 
measure is a-independent in this case, we can study both the partition function and its dual 
at ease. Let us represent the states in the fermionic Fock space as semi-infinite functions 
of \l/ r . The charge p vacuum \p) in this representation corresponds to the product: 

\p)= J] V r = y_ p+ ^_ p+l ... (5.27) 

r>— p 

To every partition k there corresponds the so-called charge p partition state: 

\ P ;k)= n = ri *- P -H*-fc 4 n * P +*-ib>= 

i=l,2,... Ki<n Ki<n 

_^ ~ ~ ~ ~ (5.28) 

JJ *_ p _i +i _ fc .|p-n> 

l<i<n 

These states form a complete basis in the space H]p\- The important fact, responsible for 
(5.25), is the variant of the boson-fermion correspondence: 

ex P ^_! b>=E§ b;k) 

^ k ^ ' (5.29) 

where factor /j,(k) is the Plancherel measure (3.12). It follows: 

a 2 

Z(a,ft,A) = A^ (0\e-x Jl A 2Lo e* J -i\0) (5.30) 
which for a = hp can also be written as 

Z(a, h, A) = (p\e-* Jl A^eR* 7 - 1 |p) (5.31) 



To make use of (5.29) we A partitions k ; , / = 1, . . . , N into a single one, K. Let 

pi, . . . ,pn be some integers. Consider the following countable set of distinct integers: 

{N{ Pl + k u -i) + l-l\l = l,...,N, ie'N} = {p + K I -I\I EN} (5.32) 
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where the sequence K\ > K 2 > . . . is defined by (5.32) and the condition that it stabilizes 
to zero. Likewise, p is determined from (5.32) and depends only on pi's: 

p = Y,Pi ( 5 - 33 ) 

Let 

ai = h(pi + pi) . 

The size of the blended partition K is expressed in terms of a;'s (see appendix A): 

1 AT 2 r, 2 n 2 

l l,i 

It is now straightforward to evaluate //g := /u(K): 

^K = ^k( a ;fr)> ( 5 - 35 ) 

and arrive at (5.25). 
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6. Af = 2 THEORY WITH ADJOINT HYPERMULTIPLET 

In a sense the most interesting M = 2 theory is the theory with the massive hy- 
permultiplet in the adjoint representation. This theory is ultraviolet finite, and is thus 
characterized by the microscopic coupling r = + ^r, and by the mass m of the hyper- 
multiplet. It is convenient to use the nodal parameter q = e 27rtT ° to count instantons. 

In this case the prepotential of the low-energy effective theory is expected to have the 
following expansion: 

.Fo(a,m,g) =7tzt af- 




The previous calculations of the coefficients fk for low values of k were done in [59]. 
One of the remarkable features of the prepotential of the low-energy effective theory of the 
theory with adjoint hypermultiplet is its relation to the elliptic Calogero-Moser system [60]. 
Indeed, in [20] an ansatz for the family of curves encoding the prepotential was proposed, 
using a version of Hitchin system [61]. This very construction of the elliptic Calogero- 
Moser system was found earlier [62] . One of the advantages of this realization is its simple 
extension to the Lie groups other then SU(N) (see, e.g. [63]). The coefficients fk of the 
prepotential of the Calogero-Moser system were calculated for k = 1, 2 in, e.g. [41]. 

6.1. Partition function 

The partition function of the theory with adjoint hypermultiplet in the O-background 
is explicitly calculable (in [9] the countour integral representation was given, the poles of 
the integral are exactly the same as those of the contour integral for the pure theory) and 
the answer is (p, = m/Tt): 




^2 [( a i - «n) 2 log (a t - a n ) - (ai - a n + m) 2 log (a L - a n + m)] + 



(6.1) 



oc 



+ $^ fc A(a, m) 




9=0 



(6.2) 




k 
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where 

Z^(a, m; h, q) = Z pert (a, m; h, q) //£(a, K) x 

TT (flj ~ Qn + m + h (j - i)) 

n Af .Aai-a n + m + h - k nJ +j- i)) (6.3) 

Z pert (a, m; K, q) = exp ^ 7^ (a/ - a n ; <?) - 7^ (a/ - a n + m; g) 
6.2. Abelian theory 

Let us consider N = 1 case first. Let n = ™. In the N = 1 case we are to sum over all 
partitions k: 

Z(a,m;R, ? ) = ? ^+A(^- 1 ) e -^(»;«) ^ g N J] ( KD £-/ ) (6.4) 

k nek ^ ^ ' ' 

The fact that (6.2) reduces to (6.4) for = 1 is a consequence of a simple identity between 
the Chern characters: 

^(ki-kj+j-i) _ e h(j-i) _ e va+h(k i -k j +j-i) _|_ e m+ft(j-i) _ e m+M(D) _ gfthO (gg) 
i<j D6k 

Another expression for (6.4) will be useful immediately: 
Z(a, m; h, q) =q£ e~^ m ^ Z^), ZQjl) = <^(m 2 -i) £ q^Z^) 

k 

Z ( ) = TT (£(k) + ^ = ? + ^ )!(£(k) + - Z ~ ^ )!(£(k) ~ Z) ' 2 x 

klw 11 (e(k)-t + fi)\(£(k)-i-fi)i(£(k) + h-i)i 2 {b - b) 

det 1 <j )J <^( k )| 



fci— fej+j— i+fi 1 
x ~ j ~ 

det l<i,j<^(k)||jZi+^ll 



Proposition 4 The analogue of the formula (5.31) for the theory with adjoint matter is: 

2 

Z(^ = q^Tr no q L °-^ V^(l) (6.7) 

where 

Z \ U.Tn I \ ^ ^ ^ 



V„(z) = : e l ^ (z) : = exp -^^^_ n £- exp f M ^,7 n £_ 1 (6.8) 

V n>0 / \ n>0 / 

<p(z) is the bosonic field (5.19), and H\ is the charge A subspace of the fermionic Fock space. 
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Proof: 

We represent the trace (6.7) as a sum over partition states: \p; k), where p = ^ is fixed. We start 
by evaluating: 

(k|v M (i)|k> = d ^<-^<^)ll^-w.^-^)H (6 . 9) 

det 1 <j ;J <^(k)||Cr-j+£(k),-j+^(k)|| 
where we employ Wick's theorem, and where for i,j 6 Z + : 

Ga =(-) 



+ sinTr^ 



- j + /J,) ixu. 

= CoefT^ (l^Y 1 = (6.10) 



1 — x / 1 — xy 



= <0|* iV^l)* i JO) 

2 2 



The second line of (6.10) follows from: 



\m - ydy + xd x ] | (i-£f ^-L- } = mi - i/)"-^ - a)""" 1 

while the last line is most easily derived from the current action on the fundamental fermions 

(Ol^y- 1 ) ^(1)^)10) = 

= expL^i(x"-y")J x(0|$(y- 1 )*(x)|0) = 

V n>0 / 

\1 — x ) I — xy 
from which (6.10) follows by employing the expansion (5.8). 

QED 

Remark: Note that an immediate consequence of the equalities (6.4) and (6.7) is following 
identity 

,k| TT Hnf - u 2 



n( i -9 B r- i =iy kl n 



h(u) 2 

n>o k nek 

(6.12) 

Z(p) = niqf- 1 

(see appendix C for notations) which interpolates between the Dyson-Macdonald formulas for the 
root systems A n . Indeed, when u. = n, the exponent on the left in (6.12) equals 

dim SU (p) = u 2 — 1 , 
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while on the right we have a sum over all partitions with no hooks of length \x, which one can easily 
identify with the sum over the weight lattice of SU(n) entering the Dyson-Macdonald formula. 



6.2.1. Higher Casimirs and Gromov-Witten theory of elliptic curve 

Just as in [10] we can consider deforming the theory by the higher Casimirs. Just as 
in [10] these are represented in the free fermion formalism by the W-generators, and the 
trace (6.7) becomes the Wi +00 character, in the presence of the "vertex" operator Vsz. 

For m = this trace becomes exactly the Wi +C<} character, as studied in [64] [65]. As 
shown in [36] this trace also coincides with the all-genus partition function of the Gromov- 
Witten theory of an elliptic curve (see also [66]). It would be nice to find a Gromov-Witten 
dual of the "vertex" operator V M , and establish precise relation of (6.12) to [36] [65] [64]. 

Note that the "flow to the pure Af = 2 theory" which is represented by the limit 
m — > oo, q — > 0, s.t. A 2 = m 2 q stays finite, leads to the dual Gromov-Witten theory of 
CP 1 with log (A 2 ) playing the role of the Kahler class. 

6.3. Nonabelian case 

Using (6.4) it is simple to express the partition function (6.2) in terms of free fermions: 

Z(a; h, m, q) = Tr^w (q Lo V^(l)) (6.13) 

f 

where 

K N) = ®i=i«IP«] 

stands for the Fock space of N free fermions, with specified charges (pi, . . . , pjv) (U(1) N 
weight, in the language of affine Lie algebras). 

We also have an expression for the dual partition function: 

Z D (Z,p;h,m,q) =Tr Hp (</ Lo V M (l) e H «) (6.14) 

Of course, the modular properties of the partition functions (6.6) (6.12) (6.14) reflect the 
S-duality of the M = 4 theory [11]. 
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6.4. Path representation 

As in the pure gauge theory case it is extremely useful to represent the partition 
function as a sum over the profiles - paths f(x). The result is: 

Z(a, m; h,q)= Z f (a, m; h, q) 

Zf(a.,m;h,q) = exp ^-^ ^ ctecfy f"(x)f"(y)-f h (x-y;q) + 

+ dxdy f"{x)f"{y) lh (x-y + m;q)+ 
+ ^og(q) J dxx 2 

(6.15) 

Interchanging x and y one sees that this is symmetric with respect tomH — m. 

6.5. Variational problem for the prepotential 

Again, just as in the case of the pure Af = 2 gauge theory we are looking for the extremum 
of a certain action functional on the space T a of profiles f(x) with fixed partial charges. 
By adding the surface tension term to the action we get a minimization problem on the 
space of all profiles, i.e. functions that satisfy 

f(x)-N\x\ = 0, |x|»0, 

and the Lipschitz condition < N. The Lipschitz condition in our case will be 

satisfied automatically, since all extrema will be convex. We also impose the total zero 
charge condition, a = 0, i.e. 

/ x f"(x) dx = (6.16) 
Jr 

In this section we set m = im. 



6.5.1. Path energy and surface tension 

Let 

T.(r\ = 

\ ar + m 



I f x 2 \ X 
L(x) = — x log ( — — ) — m arctan — (6-17) 
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We have: 



£'(*) = 5 log ( 



2 

X 



x 2 + m 2 



(6.18) 

1 1/1 1 x 



x 2 \ + im a; — im 



We introduce a generalization of the energy functional (4.28): 



£ q ,mU) =-l-f ( N + f'(*))( N ~ f'(y))L'(y - x) dx dy+ 

Jx<y (6.19) 



17VT 



x 2 f{x)dx 

R 



2 

where the first term is the leading % — > asymptotics of (6.15), and the second term 
comes from the weight (/' k '. In the limit q — > 0, m — > oo, so that A = mq^" stays finite, 
^<j,m - (this reflects the flow of the theory with adjoint hypermultiplet to the pure 
N = 2 super- Yang-Mills theory) . 



The functional which we want to maximize is: 

S q ,m{f) = -£ q ,m{f) + \ J *(/'(*)) , (6-20) 

where the last term is the usual surface tension term. Again, the functional (6.20) is 
convex, therefore every local maximizer is also a global maximizer. 

6.5.2. Variational equations 

Varying f'(x) in (6.20), and integrating by parts, we obtain the following equation 

f"{y) L(y -x)dy- a'(f'(x)) + 2nirx = const (6.21) 



~j~v=. 



which should be satisfied at any point of continuity of a'(f'(x)). The constant in (6.21) is 
the Lagrange multiplier corresponding to the constraint (6.16). At the points of disconti- 
nuity of a'(f'(x)), the equality (6.21) should be replaced by the corresponding two-sided 
inequalities. 

The function a' is piecewise constant, therefore differentiating (6.21) with respect to 
x we get the equation 

J f"(y)L'(y-x)dy = 2mr. (6.22) 
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6.6. The spectral curve 

The facets of the limit shape correspond to the singularities of the surface tension a, 
therefore the maximizer / will have N — 1 of them. This means that the support of /" 
will be N intervals (bands): 

supp/"= |J supp//' (6.23) 

l=l,...,N 

and, as in (4.57), we denote: 

supp/r = [aT,a+] (6.24) 

Consider the resolvent R(z) = Rf(z; 0, 0) introduced in (4.20). It has jump discontinuities 
along the bands with the jump equal to mf"(x). Set, by definition (cf. [67] [68]), 

G{z) = R[z-—j -R\z+ — 

This function is analytic in the domain U which is the Riemann sphere with 2N symmetric 
cuts 

U = C\\J 

i=i 

We have 

G(z) = 0(z~ 2 ) , z^oo, 
and hence we can consider the multivalued function 

G(y)dy. 



N 

im , im 
1 2 1 2 



The period of f G(y) dy around any cut is equal to the integral of the jump of G(z) along 
the cut. By construction, the jump of G(z) is ±nif"(x) and, as J a | f"(y) dy = 2 we have 

a z 

the well-defined function 



F{z)= ^n f G (y) d y mod z 



' oo 

im 



- f dy f"{y) log f - T V ) mod Z 



(6.25) 



We will show that F is an iV-fold branched cover of the standard cylinder 

C T = {vj | |S(zx7)| < -ir/2} mod Z 
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by the domain U. Clearly, F maps the infinity of U to the origin w = of C T . 



- - = 

Fig. 11 Covering finite cylinder 

The following proposition describes the behavior of the function F(z) on the both 
sides of the two cuts corresponding to any band [a~[ , a^~] 

Proposition 5 For any x 6 [a z - ,a^~], we have 

F(x+™±io}-F(x-^TW\=T. (6.26) 




Proof: 

This is, in fact, a reformulation of (6.22). Indeed, let x± denote the two arguments of F in 
(6.26). We need to show that 

+ G(y) dy = 

jl [ : \ I rr ,„ <627) 

4 « J,_ z L\v-z + ! ? v-z-'-t) 

Since the integrand is nonsingular on the domain of integration, we can change the order of 
integration, which gives 

2^ J L'(y-x)f"(y)dy = r, (6.28) 

which is precisely (6.22). 
QED 

Since f"(x) is real, we have 



R(z) = R(z) 
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from which it follows that 



G{z) = -G{z) , F(z) = F(z) . (6.29) 

Proposition (6.26) now implies that 

x± — j==Fzr/2, x G 

which means that F maps the boundary dU to the boundary dC T . It follows that F maps 
U onto C T and, moreover, since F has degree iV on the boundary, the map 

F : U -> C T 

is an iV-fold branched cover. 



Proposition (6.26) also gives the analytic continuation of the function F(x) across the cuts 
of U. The Riemann surface C of the function F is obtained by taking countably many 
copies of the domain U and identifying the point 



on one sheet with the point 



lm , - n - r - +n 

x H — — ± iu , a; G [a z , a z J 



im 

x =F zO 

2 



on the next sheet. 

By construction, F extends to the map from C to the infinite cylinder 

F : C -> C mod Z 

satisfying the property 

F _1 (c7 + r) = +im. (6.30) 

More geometrically this means that the curve C is imbedded into the total space A of affine 
bundle over the elliptic curve E T = C/Z © rZ which is obtained from C T by identifying 
its boundaries: 

C C A = {(z, vj)\z 6 C,ro6 C}/(z, w) ~ (z, vo + 1) ~ (z + m, + r) (6.31) 
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The map F is the restriction on C of the projection A — > E T . The coordinate z is the 
coordinate on the fiber of A. 

This implies, in turn, that C is the spectral curve of the elliptic Calogero-Moser system 
[60]. Indeed, the latter is the spectral curve of the Lax operator (see appendix C): 

C : Det z , n (L(w) - z) = 

L Zn = <J In p n + m— log (^iiM) + — 1 - Si n ) a ( g 7 ^— 

V 27TZ / 2tU ^ll(tJ7)^ll(9l - 9n) 

(6.32) 

and is naturally viewed as the holomorphic curve in A representing the homology class 
of N[E T ] under the identification A ~ C x E T : (z,w) = (p + mlog (6*11(^7))', ro). This 
condition, together with the behaviour near = fixes uniquely the iV-parametric family 
of curves C. The coordinates u = (u±, . . . , un) on the base of the family of curves can be 
defined using the characteristic polynomial of the Lax operator (6.32). An alternative set 
of (local) coordinates is given by the so-called action variables, which can be defined as a 
set of periods of the complexified Liouville one-differential zdw on C. In our construction 
these periods are linear combinations of a/'s, £j's and m, see below. 

Remark: For N = 1, the function F^ 1 (zu) is, up to normalization, the classical WeierstraB 
function C(^)- 

The limit shape / is reconstructed from the fact that f'(x) mod 2Z is equal to the jump 
of the function F(x) across the cut of U. 

6.6.1. Periods 

We shall now study the periods of the differential 

dS = zdw\ c = zdF(z) = ^-z (r (z - ^ - R (z + dz (6.33) 

on the Riemann surface C. 

For I = 1, . . . , N let denote the 1-cycle, which circles around the cut [aj~ ± af ± *?p]. 
Clearly: 

ivn 

dS = ai ±— (6.34) 

1 

For I = 1, . . . , N — 1, let Pi be the cycle on C joining some points x + ™ and y + *rp, where 

x E {a7,a+) , y E (aY +v a+ +l ) , 
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above both cuts on some sheet of C and below both cuts on the next sheet of C. 
By construction of C, we have 



zdF(z) = 



Pi 



x+^+iO 



f 

zdF(z) + / 

J v 



x-*f--iO 



zdF(z) = 



y-^-iO 



zdF(z) 



x+^+iO 



zdF(z) . 



x-^-iO 



By the same principle as in the proof of Proposition (6.26), we obtain 



zdF(z) = — / f"(t)dtx 



im\ , f x — t + iO\ ( %m\ , / x — t — iO 

t + — log — — + it - — log — — 

2 J \x — t — im J \ 2 J \x — t + im 



^- f (xV \t - x) + L(t - x)) f" \t) dt 
2m J r 



(6.35) 



= — 6+i + const 
2m 

thanks to (6.21)(6.28) and the fact that a'(f'(x)) = 6 because x lies on the /'th band. We 
summarize: 





- 6+1 = 2ttz 


I dS 

hi 






dS 


a L 








.F(a 


m, g) = -£ q 


,m(/*) 



(6.36) 
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7. ON ANOTHER MATTER 

The results of the previous chapters can be generalized in various ways. One can 
incorporate hypermultiplets in arbitrary representations. One can add a tower of Kaluza- 
Klein states. One can study quiver gauge theories. We shall only sketch the results for 
the fundamental matter, or for five dimensional theories. The other cases will be treated 
in future publications. 

7.1. Hypermultiplets in the fundamental representation 

In this section m denotes the vector of masses: 

m = diag(mi, . . . ,m k ) 

We will denote by f = 1, . . . , k the flavour index. 

The partition function of the theory with k matter hypermultiplets in the fundamental 
representation of U(N) is given by the following sum over colored partitions [9]: 

Z(a, m; % A) = A( 2N ~ k W Z c (a, m; h, A) 
k 

Z £ (a, m; h, A) = Z*"*(a, m; h, A) J[ L JL_ ^ _ ^ ^( a , 
Z pert (a, m; ft, A) =exp K] 7 ^ (a/ - a n ; A) + 1% (a t + m f ; A) 

V l,n l,f 

As before, the partition function (7.1) can be written as a sum over profiles / = f a £ as 
follows 

Z(a, m; ft, A) = ^ Z/(a, m; ft, A) 
/ 

Z/(a, m; ft, A) = exp j dxd V f"{ x )f"{v) In (x - y; A)+ ( 7 . 2 ) 

f dx f"( x ) ^ h ( x + mf; A ) 

Observe that the action in (7.2) is the old action (4.25) evaluated at a modified function 
/ such that 

/" = /"-£*(* + "*) (7-3) 
f 
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Therefore, the variational problem for the limit shape / is essentially the variational prob- 
lem solved in Section 4, the only modification being that the maximizer /* is required to 
have corners at the points x = —nif. 

We will solve this problem under the assumption that rrif are real and \rrif\ ^> 0. 
Without loss of generality, we can assume that the number k of flavors is exactly 2N — 1 
since letting some rrif go to infinity reduces the number of labeled. 

Let A be the half-strip 

{\^(w) \ < N,%w) > 0}, 

but now with two kinds of vertical slits. As before, we have N — 1 slits going from the 
points 

-N + 21 + irn, l = l,...,N-l 
down to the real axis. Additionally, we introduce 2N — 1 slits from the points 

-N + l + irfi, I = 1,...,2N-1 

going up to the infinity. Let $ be the conformal map from the upper-half plane to the 
domain A as in the following figure 



m 3 m 2 m, oo m 5 m 4 




Fig. 12 The map $ for iV = 3 

The numbers r\\ are, as before, the parameters of this map, as well as the number A. The 
numbers rji are uniquely fixed by the requirement that 

$ -1 (oo) = {mi, . . . , TO2JV-1, oo} • 
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Let P(z) be a monic degree iV polynomial with all real roots. Suppose that A is 
sufficiently small and the mf's are large enough. If 



mi < . . . < m r < < rn r+ i < . . . < m 2 7v-i 

then, as before, one checks that the map 

2 w 
= — In —= - N + r + 1 , 

where w is the smaller root of the equation 

A Pi ) 27V-1 

» + " = ^, Q(,)= I[(, + m f ), 



(7.4) 



(7.5) 



is the required conformal map. Its construction is illustrated in the following figure which 
is a modification of Fig. 7: 






2A 1/2 





2A m 



Fig. 13 The map $ as a composition of maps for N = 2 
The branch of the logarithm in (7.4) is chosen so that 

Slog (w) — » , z — » oo . 
We now claim that the maximizer /* is given by the same formula as before 

A(x)' = 3?$(x), (7.6) 
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which simply means that the nontrivial part of the maximizer /* is determined by the 
formula 

/*(x)' = 9fc$(x), |x| < min{|m f |}, (7.7) 

whereas f*(x) = N\x\ outside of this interval. 

Indeed, f*{x)' clearly has the required jumps at the points x = —nif. The equation 

(X/*)' = -7rS$ 

is checked in the exact same way as before. In particular, the asymptotics 

<&(z) ~ — In — + real constant , z — > oo (7.8) 

TVl Z 

agrees with the fact that from (4.61) we have 

(x£)'(z) =lnA-ln|z| +0 Q , z -> oo . (7.9) 

This verifies the predictions of [3] [4] . If 2 N = k then the curve will be covering elliptic 
curve, and again can be explicitly constructed. 

7.2. Five dimensional theory on a circle 

Consider five dimensional pure super symmetric gauge theory with eight supercharges, 
compactified on the circle S 1 of circumference (3. In addition, put the fi- twist in the 
noncompact four dimensions. The resulting theory is a deformation of (2.1) by /9-dependent 
terms. 

The theory was analyzed in [14] [9] and the result of the calculation of the partition 
function is the following: 

Z(a;^,/3,A) = ^(/3A) 2Ar l £ IZ £ (a ; n,/3) 
k 

%(a;M) = ^ eH (a;M) ^(a,/3,ft) 

2/ a k\ - TT sinn f( a z ~~ a n + H k i,i ~ k n ,j + j - i)) (7.10) 

liAa,p,n)— II . . a, 7TT- ^ 

(i,0*(»,j) sinhf(a ; -a n + ft(j-z)) 

^ ert (a; h, (3) = A^exp £ 7R (a, - a n |/3; A) 

l,n 

Note that we included a power of (3 in A in (7.10) in order to have a simple four dimensional 
limit (3 — > 0, A finite. This "renormalization" group relation was actually derived in [6] [14]. 
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7.2.1. Path representation 

The partition function (7.10) is easily written as a path sum: 



Z(a;M,A)= E ex pf-l/ f"(x)f"(y) lh (x-y\(3;A)) (7.11) 
7.2.2. Prepotential 

From (7.11) we derive, as before, the equations on the critical path, which, in turn, 
determine the prepotential. Assuming (3 > 0, we have 

^(a.jS.A) = jZpif*) , (7-12) 
where /* is the maximizer and the functional £p is defined by: 

W) = / < (^ + f(x))(^-/'(y))log^ sinh ^_^^ . (7.13) 

Again, we modify the functional (7.12) by a surface energy term to obtain the following 
total action functional: 

S{f) = \SftU) + \J <T<J')dx (7.14) 

We now introduce a five dimensional analogue of the transform X/ (see appendix A for 
the definition of the function 70 (x; 13)): 

t x /3/] ( x ) = \i dy sgn(y-x)y (\y-x\;p)f"(y). (7.15) 



Note that the relation 



implies that 



7o(x;/3)'' = log^sinh^, (7.16) 



[X,/]' (x) = -\j dy log (jL sinh ^_A^ f> {y) . (7.17) 

The equations on the minimizer for (7.14) are formulated exactly as in the Proposition 
(4.45) with the replacement of X/ by X^/. 

The solutions to these equations can be constructed as follows. Let $>(z) be the 
conformal map from the horizontal strip 

< %z < ^ (7.18) 
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to the slit half-strip A with N — 1 slits of the form (4.48) and an additional vertical 
semi-infinite slit 

JJtu = , > rj , 

as in Section 7.1 above. We normalize $ by the requiring that it maps the two ends of 
the horizontal strip to the two ends of A. This fixes it up to precomposing with an overall 
shift by a real constant. 
The exponential map 

z X = e pz 

maps the strip (7.18) to the upper half- plane, mapping R to R>o- Therefore, we can use 
the results of the previous section to conclude that $ has the form 

2 

$(z) = N -\ logic? , 

TTl 

where w is the smaller root of the equation 

X~%P N (X) = ((3A) N (w + , X = e Pz . (7.19) 

Here P/v is a monic polynomial of degree N and logu> is normalized, as usual, by the 
requirement that 3 log w — » as z — » +oo along the real axis. 

In the form (7.19), the parameter rj of the map $ becomes a function of A and other 
parameters. The equation (7.19) defines the spectral curve of the relativistic Toda chain, 
as predicted in [14]. 

Applying the Schwarz reflection principle to any vertical part of the boundary of A 
and taking the resulting function modulo 1, we obtain a 2 N- fold map from the horizontal 
cylinder of circumference 2tt / (3 with cuts along the real axis to the unit half-cylinder, 
as shown in the following figure 




Fig. 14 The map $ mod 1 
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Here each cut covers 2-to-l the finite boundary of the half-cylinder, while both ends 
of the cylinder cover iV-to-1 each the infinite end of the half-cylinder. 

As before, it follows that &'(z) extends to a 27ri//3-periodic function in the entire 
complex plane with N cuts along the real axis, periodically repeated. The values of $>'(z) 
on the two sides of the cut are real and opposite in sign. We have 

/ Nfix \ 2N 
<S>(x + iO) = ± [i— — + N )+— log/3A, x -> ±00 , (7.20) 

\ TV I 1X1 



whence 



N3 

3>'(x + iO) = =F — t- , x^±oo, (7.21) 

TVl 



It follows that the function <&'(z) satisfies the following integral equation 

&( z ) = — [ m'(x) coth (p ( X ~ Z A dx. (7.22) 
Moreover, by comparing the asymptotics, we can integrate (7.22) to 

§( z ) = N+ - [ m'(x) log ( 4t sinh ^ z ~ X A dx . (7.23) 
We claim that the maximizer f+ is given in terms of the map $ by the usual formula 

fl(x) = m(x + io) . 

Indeed, from (7.23) and (7.17) we conclude that 

(XM' (x) = ~^(x) . (7.24) 

The difference with the previously considered cases lies in the fact that to reconstruct the 
values of £/'s (and a/'s) we have to integrate the multi- valued differential 

dS 5 = ± log (X)^- 
p w 

on the curve (7.19), as indeed proposed in [14]. Also note that domain where the map $ 
is defined is different. This leads to the new form of the Seiberg-Witten differential. 
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7.2.3. Free field representation 

Again we start with N = 1 case. Let Q = e^ h . In this case we are to sum over all 
partitions: 

/ \ 2 

z(a,h,p,A) = j2(-^ l m) m n 



□6k 



2sinh(«a) 



oo 



(7.25) 



exp^ 



{(3K 



,2n 



jf=! An sinh 2 (^) 



The free field representation of the sum (7.25) is most simply done using the vertex oper- 
ators: 

r ± (/3,ft)=exp £ ndlnn^ (7-26) 

±n>0 ^ ' 

Then 

Z(a,h,P,A) = (0|r + (/3,n)(/3A) 2L °r_(/3,^)|0) (7.27) 
and, for A > 1 we would get: 

Z D (Z;p;P,h,A) = (p\T + (P,h)(PA) 2L "e u ^.(P,fi)\p} (7.28) 
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8. CONCLUSIONS AND FUTURE DIRECTIONS 



8.1. Future directions and related work 

In this paper we extensively analyzed the partition functions of various H = 2 super- 
symmetric gauge theories, subject to the O-background of M = 2 supergravity. In all cases 
the partition function was identified with the statistical sum of an ensemble of random 
partitions with various Boltzmann weights. The grand canonical ensemble corresponds to 
the fixed theta angle in the gauge theory. In the thermodynamic limit the statistical sums 
are dominated by a saddle point, a master partition, from which one can extract easily 
various important characteristics of the low-energy effective gauge theory, such as the pre- 
potential of the effective action, which is identified with the free energy per unit volume 
in the dual statistical model, provided we relate the volume to the parameters ei, e 2 of the 
O-background, via 



ei£2 

The particle number is clearly the instanton charge. So the gauge theory partition function 
corresponds to the grand canonical ensemble, with log (A) being the chemical potential. 
The average particle number per unit volume, the density, is easy to determine: 

e i e 2 d ^ inst 




2N Z inst (a, ei, e 2 , A) <91og (A) 

(recall (4.5)). In the region the right hand side of (8.1) is clearly very small, it is of 
the order: 

A 2N 



2 „ 9 



' a 2(7V-l) 

where M\y is the typical VF-boson mass, and g e ff(Mw) is the effective coupling at this 
scale. The dependence of the instanton density on the coupling is the typical one in the 
dilute gas approximation (cf. [23] [69]). In the semi-classical region there are very few 
instantons per unit volume. 

Of course, a cautious reader may wonder about the dimensionality of the instanton 
density (which normally should be (mass) 4 ). The point is that the noncommutative regu- 
larization which we used in the calculations introduces a scale ~ v 7 ©, which is responsible 
for the dimensional transmutation we see in (8.1). Indeed, the physical volume, occupied 
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by k instantons, sitting on top of each other, is roughly k<S) 2 . Such instanton "foam" [70], 
or perhaps, liquid, may well be related to more complicated phenomenological pictures, 
say instanton liquid [71]. Moreover, as we go deeper into the moduli space of vacua, the 
diluteness of the instanton gas also must cease to hold. However, the analytic properties 
of the partition function (3.2) are powerful enough to uniquely fix it by the instanton gas 
expansion in Uoo . 

Note, that one of the interpretations of our result is the existence of the master field in 
the M = 2 gauge theories, which, however, is working only for special (chiral) observables, 
unlike the master field of the large N gauge theories [72]. The latter, however, is much 
more elusive. 

In our story, the master field is constructed as follows. First of all, the labeled parti- 

— * 

tions k we sum over are in one-to-one correspondence with the (noncommutative) gauge 
fields, describing certain arrangements of the nearly point-like instantons, sitting nearly on 
top of each other. The "master" partition with the profile /* corresponds to statistically 
most favorable configuration. 

We extensively exploited in this paper the fact that O-background acts as a box, 
similarly, in some respects, to the AdS space for the supergravity theory. One may wonder, 
why couldn't we use the more traditional ways of regularizing infrared divergences of the 
gauge theory. In fact, the number of options is rather limited. One may study gauge 
theories on the compact four manifolds. To preserve supersymmetry one needs to turn on 
certain non-minimal couplings, which effectively twist the gauge theory. However, in these 
approaches one cannot learn directly the properties of the gauge theory in the infinite 
volume, as all vacua are averaged over, and either most of them do not contribute to 
the correlation functions of the chiral operators (which is the case for the simple type 
manifolds), or (for the manifolds with < 1) the contribution of the various vacua is 
related to the prepotential in a complicated fashion [73] [74] . It is out of such attempts to 
regulate the Af = 2 theory that one naturally arrives at the concept of the O-background. 

The topic which we completely neglected in our paper is the application of our for- 
malism to the theories with Af = 1 supersymmetry. Recently, there was a lot of excitement 
related to the exact calculations of the effective superpotentials of the theories, obtained 
from Af = 2 theories by the superpotential deformations [8] [75] [76] . It should be straight- 
forward to apply our techniques in these setups as well. One encouraging feature of our 
formalism in the case of pure Af = 2 and softly broken Af = 4 theories, is the striking 
similarity of the expression we have for the partition functions and the matrix models 
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calculating the effective superpotentials of the Af = 1 gauge theories with an additional 
adjoint chiral multiplet, or with three adjoint chiral multiplets, with specific superpotential. 

The similarity is the measure — in the first case it is the regularized Vandermonde 
determinant of the infinite matrix $ with the eigenvalues ai + h(ku — i). In the second 
case it is the ratio of the determinants: 

Det'(ad8) 
Det(a<i$ + m) 

The fact that the eigenvalues concentrate (for small %) around a/'s also has a matrix model 
counterpart. Namely, to get the theory with U(1) N gauge symmetry one adds to the U(N) 
theory a tree-level superpotential with extrema near ai, . . . , o/v- In the limit of vanishing 
superpotential one is left essentially with the Haar measure on matrices, together with the 
prescription to distribute eigenvalues near the extrema of the phantom of the superpotential 
[8]. 

The differences are also easy to see: on our case the matrix is infinite, but the eigenval- 
ues are discrete and are summed over. In the Dijkgraaf-Vafa case, the matrix is finite, but 
the eigenvalues are continuous and integrated over. We have % — > 0, they have N — > oo. 
The issue is tantalizing and is under investigation [77]. 
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Oh, East is East, and West is West, and never the twain shall meet, . . . 

Another issue which we didn't touch much upon is the string theory dual of our calculation. 
There are several points one may want to stress. 

First of all, the summation over partitions that we encountered is very similar to 
the summation over partitions one encounters in the two dimensional U (A) Yang-Mills 
theory, in the large A limit. In particular, when working on a two-sphere, one finds a 
master partition [78], and moreover, for the values of the 't Hooft coupling constant the 
profile of the master partition has the facets, just like ours. 

Of course, technically speaking, the two dimensional YM theory corresponds to the 
U(l) theory in our case, and the 't Hooft coupling there corresponds to the higher Casimir 
coupling T3 in our game. But the phenomena have similar origin. Now, the two dimen- 
sional Yang-Mills theory has a dual closed string representation [79] , and so does our four 
dimensional partition function. And again, just like 1/A played the role of the string 
coupling constant, % plays this role in the four dimensional story. 

In the U(l) case all this is more or less well-known [36] [10] by now. However, the full 
string dual of the U (A) theory is yet to be discovered (see [10] for the discussion, see also 
[80] [81]). 

We should add that from the point of view of Gromov-Witten theorists the dual type A 
topological string theory, whatever it is, must be rather simple. Most of them should have 
one dimensional target spaces. The rough correspondence between the gauge theory and 
the dual GW theory states that the U(N) theory with g adjoint hypermultiplets (which 
is not asymptotically free/conformal for g > 1) is dual to the GW theory of A copies of 
the Riemann surface of genus g. The words "A copies" still do not quite have a formal 
meaning. 

In the large A description of the two dimensional Yang-Mills theory a prominent role 
was played by the formalism of free fermions [53]. Of course, the chiral fermions of [53] 
are ours \1/ (somehow four dimensional gauge theories do not see the anti-chiral sector). 

One of the exciting problems is to understand better the nature of these chiral fields 
in the string realizations of M = 2 gauge theories. The conjecture of [9] was that they 
arise as the modes of the chiral two form propagating on the worldvolume of the NS5/M5 
brane, trapped by the Q-background. In the bosonized form, these can also be mapped to 
the truncated version of the Kodaira-Spencer field, propagating along the Seiberg-Witten 
curve. We are planning to investigate these issues in the future. 
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Yet another exciting area of research, revolving around our partition functions is their 
two dimensional (anyonic) interpretation in the v 7^ 1 case. The connections to the Jack 
polynomials, which are the eigenfunctions of the Sutherland many-body Hamiltonians with 
v{y — 1) coupling constant, suggest strongly some relations to the physics of the quantum 
Hall effect, and also to the theory of analytic maps [82]. The latter also seem to be 
responsible for the higher Casimir deformation (4.73). 

Seiberg-Witten curve itself seems to arise as some quasiclassical object in the theory 
of some many-body (N — > 00) system. Indeed, the space z, w where it is embedded, can 
be viewed as a one-particle phase space. The particles are basically free for v = 1, and 
have some phantom interaction for v 7^ 0, 1, which only reveals itself in the generalized 
Pauli exclusion principle (Haldane exclusion principle). They fill some sort of Fermi sea, 
bounded by the Seiberg-Witten curve. Perhaps, the relation to the random dimers [47] 
will also prove useful in understanding these issues. Another possibly useful feature of 
our energy functional is its seimple relation to the two dimensional local theory of a free 
boson. Indeed, the bi-local part of (4.36) is nothing but the induced boundary action, 
with p(x) being the boundary condition in the theory of the free boson, living on the 
upper half-plane. Of course, when the subleading in h corrections are taken into account 
the theory will cease to be free. It is also interesting to point out that our extremizing 
configuration f+(x) corresponds to some sort of multiple D-brane boundary state, with N 
D-branes located at x ~ a/, I = 1, . . . , N. 

It is of course very tempting to develop these pictures further, connect instanton 
gas/liquid/cristall to the (Luttinger?) liquid of the dual anyons, and learn more about the 
properties of the fivebranes from all this. Perhaps, the results of [83] [84] will prove useful 
along this route. 

Finally, the theory with adjoint matter should be closely related to conformal field 
theory on an elliptic curve. We have uncovered some of the relation in the partition 
function being interpreted as a trace in the representation of a current algebra. We should 
also note that the elliptic Calogero-Moser system, whose spectral curve, as we showed 
above, encodes the quasiclassical/thermo dynamic limit of the partition function, allows a 
nonstationary generalization, related to KZB equations, which, conjecturally, governs our 
full partition function [85] [86] [87] [88] [89] [90] . 
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8.2. Summary of the results 

In this paper we have advanced in the study of the vacuum structure of gauge theories. 
We considered M = 2 supersymmetric gauge theories in four dimensions, and have subject 
them to the so-called O-background. In this background one can calculate exactly the 
partition function, in any instanton sector. The result has the form of a statistical sum 
of the ensemble of random partitions. The limit where the O-background approaches flat 
space, is the most interesting for the applications, as there one is supposed to learn about 
the properties of the supersymmetric gauge theory in flat spacetime. In terms of the 
random partitions this is a thermodynamic, or quasiclassical limit. The flat space limit of 
the free energy coincides with the prepotential of the low-energy effective theory. 

We have evaluated the sum over instantons, by applying the saddle point method, and 
thus have succeeded (to our knowledge, for the first time in the literature) in producing 
the all-instanton direct calculation of the prepotential. 

In all cases considered: pure gauge theory, theory with the matter hypermultiplets, 
five dimensional theory compactified on a circle - the saddle point corresponds to some 
master partition, which is the analogue of the eigenvalue distribution in the theory of 
random matrices. And in all cases one can encode the solution in some family of algebraic 
curves, endowed with a meromorphic (sometimes multivalued) differential, whose periods 
contain the information about the prepotential. 

We have also found interesting representation of the full partition function (which, 
in addition to the prepotential, contains also certain higher gravitational couplings T g of 
the gauge theory) as a partition function of the theory of chiral fermions/bosons on a 
sphere (for the pure gauge theory), on a torus (for the theory with adjoint matter), or 
some (/-analogue thereof (for the five dimensional theory). 

We have also uncovered numerous puzzling relations between various seemingly unre- 
lated topics which leave a lot of work for the future. 

Disclaimer. Opinions presented in this paper do not necessarily reflect author's point of 
view. There are about ~ 2%f misprints in this paper. We are however confident that the 
most of them cancel each other. 
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Appendix A. The function 7^ (x; A) 
A. 0.1. Free case: e 2 = — e\ = h 

The function 7^ (x; A) is characterized by the following properties: 

1. Asymptotic expansion for % — > 0: 

00 

7,(x;A) = 5> 2s - 2 7 s (x) (A.l) 
3=0 

2. Finite-difference equation: 

ln (x + h;A)+ lh (x - h; A) - 2 7 ^ (x; A) = log (A.2) 

The conditions (A.1)(A.2) specify 7^, (x; A) uniquely up to a linear function in x. All 
the terms r y g (x), g > are uniquely determined: 

3 2 

"4* 



7o(x) = ^x 2 log 
71(X) = -T2 l0g (l) 



_B 1 

7g(x) = 2^ 2 -2)^ - ^ 



where -B n 's are the usual Bernoulli numbers: 



- 1 " 2s n\ 

n=0 



The function 7^ (x; A) is closely related to the gamma function: 

7* (* + §; a) - 7 » (* - §; a) = log ( 1 «f r (i + 1 ) ) (A.4) 



Another definition of the function 7^ (x; A) is through the zeta-regularization: 



7»(*;A) = - 



A s A 00 eft s e~ te 

r( s )7 * * ( e »*-i)(e-»*-i) ( - 5) 

s=0 
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Remark: Note, that 

lh (0;A) = -^ 



The function 7a (x; A) for A = 1 arises as the free energy of the c = 1 string with h 
being the string coupling, and x the cosmological constant. It is also related to the free 
energy of the topological type A string on the conifold (see also below). The adepts of 
the applications of matrix models in susy gauge theories [8] praise yet another property of 

In (x; A): 

log (VoH7(A0) = 71 (AT; 1) (A.6) 
A. 0.2. Anyon case: general e 1? e 2 



7 £1 ,e 2 (*;A) = ^ 



r(s) J t (e £ i* - l)(e £ 2* - 1) 

s=0 



In the case e 2 = — ei = 7i this reduces to 7^, (x; A): 

7^ (x; A) = 7_^ (x; A) 

The main difference equation: 

76i,6 2 (>; A) + 76i,62 (x - ei - e 2 ; A) - 7 £l)£2 (x - ei; A)-7 ei , e2 (x - e 2 ; A) = 

A\ (A.8) 



x 



log 
Reflection: 

2x 4- 

7 _ ei , e2 (x;A) = - 7eije2 (x-ei;A) __log(A) (A.9) 

For v e Q the function 7 eii e 2 (x; A) can be related to 7^ (x; A). Suppose 

ei = --, e 2 = -, N. (A.10) 

V Q 

Then 

p ~ lq ~ l ( % \ 

76i,62 (z; A) = ^ T» ( x + — - ?») ; A J (A- 11 ) 

i=o j=o \ PQ / 
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7. 



ei ,£2 



A. 0.3. Trigonometric analogue of 7 €lj£2 (a;; A) 

The natural generalization of (A. 7) is the function: 

x 7 n=l 

(A.12) 

which obeys: 

The main Q-difference equation: 

7ei,e 2 (^1/3; A) +7e!,e 2 (a - ei -e 2 |/3;A) -7 ei , ea (x-ei|/3; A) -7 ei>ea (x-e 2 |/3; A) = 

" log (|v sinh T 

(A.13) 

We shall only use in this paper the special case e 2 = = Ti, where we get the 
function 

[3x 3 x 2 . fix ^1 e~P nx 



7 „(x;ftA)=^ * log(/3A )_|£ + E I 



12ft 2 2ft 2 24 ^ n (e-^ - l)(e^ - 1) 

n=l 

oo 

= ^ 7fl (^;/3)^- 2 

7o(x;/3) = -ylog(/?A) + - -1 Li 3 (e"^) (A ' 14) 
7i(x;/3) = -y^log ( 2sinh— J 

Note that up to the terms of instanton degree zero the function 7^ (x\fi; A) coincides 
with the all-genus free energy of the type A topological string on the resolved conifold, 
with fix being the Kahler class of the P 1 , and /3ft the string coupling [91]. Another 
Gromov-Witten interpretation is via the local F 1? with the Kahler class of the base 
P 1 being log (/?A) (considered to be big), and the fiber P 1 with the Kahler class fix 
[6]. 
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Appendix B. Partitions, charges, colours 

Partition: a nonincreasing sequence of nonnegative integers, stabilizing at some point 
at zero: 

k : ki > k 2 > ■ ■ ■ > k n > = k n+1 = k n+2 = ... 
n = £(k) is called the length of k, 

|k| = Ylii i s caue d the size of the partition, k^s are called the parts of the partition. 




Fig. 15 Box, hook, .... 



The parts of a partition are labeled by i,j = 1,2, ... , 
For the partition k: 

□ = (i,j) e k 1 < i, 1 < j < ki 

Dual partition: k, 

(i,j) ek^(j,i) e k 

i.e. ki = < kj} 

Hook-length hij = h(n) of the box in the Young diagram of the partition k: 

hi j = kj + ki % j ' + 1 

Coloured partition: k, the A^-tuple of partitions: 

k= (ki,...,kjv), 

individual partitions are denoted as: 

k ; = (k tjl > k ij2 > ■■■> h jni > ki, ni+1 = = ...), 
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|k| = J^feji 

l,i 

Charged partition: (p; k) , the set of nonincreasing integers k\ = hi + p, where 

k = (fci > k 2 > . . .) 

is a partition, and p G Z. The limit = p is called the charge. 
Blending of coloured partitions: given a vector p = (pi, . . . ,pn), with 

I 

and an iV-tuple of partitions k, we define the blended partition K, as follows: 

{Ki - i | ieN} = {W(fc K -i + pi) + J-l | I = 1, . . . ,N, zGN} (B.l) 
B.0.4. Power-sums 

To the charged partition (p; k) it is useful to associate the following shifted-symmetric 
generating function, analytic in t, with a single pole at t = 0, defined for dit > by the 
series: 

oo 

P P;k (t) = ^e^ +fc - l+ 2) (B.2) 

i=l 

The expansion of p p; k(£) near t = contains information about the charge, the size of k, 
etc: 

P P ;k(*) = \ + V + 1 + |k| - ^ + . . . (B.3) 

• 

Given N partitions k/, / = 1, . . . , AT, and the charges G Z, we associate to them the 
generating function: 

Ptf M = E e'( K(K+ ""-" + '^) = JJ e«« p pl ,,(M) (B.4) 

Z,i I 

which corresponds to the blended partition K of charge: 

P = ^2pi = ^2pi (B.5) 
l i 
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and size: 

N 2 - 1 

iki = j> 7 = ny: m + M - \v 2 - \^ (b.6) 
i i 

where 

Pl=Pl + Pi 



The function f p -k(x) is related to p p; k(*) by the integral transformation: 

/;. k (x) = -- / dt e"'** sin ^ p p;k (ztft) (B.7) 
7TZ j R z 

The resolvent (4.20) is related to p p; k(£) by another integral transformation: 

i r°° th 

R(z\e u e 2 ) = : / dt e~ tz sinh — Pp;k (tft) (B.8) 
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Appendix C. Theta- function 

For completeness we list here the relevant formulae for the odd theta function, which 
we use in (6.32). 



From these formulae one easily concludes that the Lax operator L(w) (6.32) is the mero- 
morphic Higgs operator in the rank N vector bundle over the elliptic curve E T , twisted by 
the one dimensional affine bundle, which makes its spectrum to live in the affine bundle 
as well. 

In (6.12) we use Dedekind eta- function: 




6>ii(ro + l;r) 
6>ii(ot + t;t) 



6h(zu;t) 

•e- 7ri ( 2w+T ^n(ti7;r) 



(C.l) 



oo 



r}(q) = - q n ), q 



(C.2) 
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